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EARTHQUAKE RESPONSE OF FIXED OFFSHORE STRUCTURES 



Owen Major Kirkley, Ph.D. 

Department of Civil Engineering 
University of Illinois at Urbana-Champaign , 1973 

The response of fixed offshore structures, such as those used for 
oil drilling and production purposes, to earthquake excitation is investigated. 
An approximate method of computing maximum response quantities in the structure 
is developed. The approximate method uses superposition of classical normal 
mode responses as determined by response spectrum techniques. The structure 
is modeled as a lumped parameter system, the earthquake is represented by the 
deterministic time history record of the ground acceleration, and the hydro- 
dynamic forces are considered to be those computed from Mori son's equation. 

A single degree of freedom oscillator with nonlinear drag character- 
istics is subjected to various representations of ground motion, including 
earthquake motions. The existence of a high frequency tail in the response 
spectrum is shown, and a method for linearizing the nonlinear hydrodynamic 
drag force is developed and shown to be satisfactory. A method is developed 
for conservatively estimating the magnitude of equivalent linear hydraulic 
damping in the system based on velocity characteristics of the ground motion. 
Rules are developed for the construction of idealized response spectra for 
several different response quantities. The effects of the nonlinear hydro- 
dynamic drag force on the inelastic response of the single degree of freedom 
oscillator are investigated. 

The principles and approximate methods developed during the study 
of the single degree of freedom system are applied to the multi -degree of 



freedom model of the actual structure. The linearization technique is ex- 
tended to the multi -degree of freedom system in such a manner as to insure 
the existence of classical normal modes in the equations of motion. A series 
of multi -degree of freedom systems is subjected to one of the earthquake 
records and the linearization technique is shown to be satisfactory in com- 
puting the response quantities of interest. The response of a series of 
offshore towers to an earthquake record is computed by the superposition of 
modal responses as determined by both exact spectra and idealized spectra. 

The approximate method considers the effects of the structure-fluid inter- 
action by conservatively estimating the amount of linear hydraulic damping 
developed as a result of this interaction. 
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CHAPTER 1 
INTRODUCTION 



1.1 General 

The design of fixed offshore structures is frequently based on the 
static forces resulting from the loads imposed on the structure by wind, 
waves, ocean currents, and earthquakes. Design on the basis of static forces 
is satisfactory when the response of the structure is not magnified by the 
dynamic nature of the loads. However, the characteristics of the structure 
and the forces acting upon it are often such that the response of the struc- 
ture to wind, waves, currents, and earthquakes is dynamic and magnified greatly 
over the static response. In these cases, design should properly be based on 
the dynamic response of the structure. Various methods are available for cal- 
culating the dynamic response of a land-based structure to earthquake forces, 
but similar calculations for offshore structures are complicated by the pres- 
ence of nonlinear hydrodynamic forces acting on submerged portions of the 
structure. Further, these nonlinear forces are dependent not only on the 
physical characteristics of the structure, but also on the actual response of 
the structure to the forces acting on it. 

Many offshore structures are currently in use, with the most notable 
use being that of oil drilling and production platforms. While platforms in 
water depths of up to 300 feet are fairly common now, the search for new 
sources of oil and gas may result in many offshore structures situated in 
greater water depths. The location of oil resources off the coast of California, 
in the Gulf of Alaska, and other portions of the circumpacific earthquake belt 
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as well as other regions of seismic activity emphasizes the need to be able 
to predict the response of offshore structures to earthquake motions. Drill- 
ing and production platforms represent a very large investment in cost, and 
large numbers of personnel are involved in the operation of such a structure. 
Failure of one of these platforms could result in many lives being lost as 
well as the loss of investment both from the cost of the platform and the loss 
of production time. 

1.2 Description of the Problem 

An offshore structure of the type mentioned in the preceding sec- 
tion is shown schematically in Fig. 1.1(a) and is a continuous system with 
an infinite number of degrees of freedom. Such structures are space frames 
composed of circular pipe sections generally welded together, with diagonal 
bracing as shown in the figure. The top portion of the structure may consist 
of more than one level and contains the heavy equipment necessary for oil 
drilling and production and often contains facilities for berthing and messing 
of the crew. The operating and production portion of the structure is gener- 
ally quite heavy. The research reported in this thesis is based on a lumped 
parameter model 1 of the structure with mass, stiffness, forces and other 
characteristics of the structure lumped into a discrete number of nodes, as 
shown in Fig. 1 .1 (b). 

The equations of motion for the system shown in Fig. 1.1(b), sub- 
jected to ground motion in quiescent water, are, in matrix form: 

[" m -]{x s } + C c s ]{x s - x g } + [K]{x s - Xg } = {P} (1.1) 

★ 

Numbers refer to List of References at end of text. 
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where 


FmJ 


= 


mass matrix of the system 




[C $ ] 


= 


structural damping matrix of the system 




[K] 


= 


stiffness matrix of the system 




% 


= 


absolute acceleration of node j 




• 


= 


absolute velocity of node j 






= 


absolute displacement of node j 




• 

X 3 


= 


ground velocity 




x g 


= 


ground displacement 



The symbol denotes a matrix in which the off-diagonal elements are zero. 

No structural dampers are shown on the model, but the damping matrix [C s ] can 
take a variety of forms as will be discussed in a later chapter. 

The hydrodynamic forces acting at each node of the system are con- 
sidered to be those computed from Morison's equation 2 for circular cylinders 
as: 

{P} = - IT* m J{y - r K d -J{x s |x $ | > (1.2) 

where 

['m -J = added mass matrix 3 = pC^C^VJ 
p = mass density of the fluid surrounding the structure 

C y = empirical inertia coefficient at node j 

j 

V. = enclosed volume at node j 
J 

["KpJ = drag factor matrix = •^•pIT'CpAJ 
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Cp = empirical drag coefficient at node j 
J 

A. = projected area at node j 
J 

If the coordinate transformation, {Z} = {x s - x g }, is made, where 
{ Z} is a vector of a motion quantity of a node relative to that of the ground, 
Eqs. 1.1 and 1.2 can be combined to yield: 

FM J{Z} + [C s ]{Z> + [K]{Z> = -['Mj{x g > - [~ Kp J{x $ | x $ | } (1.3) 

where 

l 

['M J = virtual (or total) mass matrix = C'mJ + ['mJ 
x = ground acceleration 

Equation 1.3 differs from the equations of motion for a land-based 
structure subjected to earthquake forces in two important aspects. First, 
the presence of the added mass term causes the structure to possess natural 
frequencies and mode shapes which are different from those of a structure with 
the same stiffness and only the structural mass ["m J. Secondly, the hydro- 
dynamic force term [~ Kp-J{x s |x s |} is not present in the conventional land-based 
earthquake problem, and this nonlinear hydrodynamic force causes the major dif 
ficulty in solving Eq. 1.3. 

If the nonlinear hydrodynamic force ["Kp ~J{x s |x s |} can be converted 
into an "equivalent" linear hydrodynamic force, say [~Kpp-J{x s }, Eq. 1.3 can 
be written as 

* 

D-MOtZ} + [C S ]{Z> + [K]{Z> = -C'MOfy - r-K DL J{x s } (1.4) 

Since {Z} = {x - x }, Eq. 1.4 can be rewritten as 
s 9 
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FMjm + [C S ]{Z} + C-K dl J{Z} + [K]{Z} 

= -['Mj{x g } - rK D[ _J{x g } (1.5) 

Equation 1.5 differs from the equations of motions for a land-based multi- 
degree of freedom damped structure of mass ["M-J, subjected to a base exci- 
tation, in the linear drag force K|^-J{x g }. 

One well-known and convenient method of solving simultaneous linear 
differential equations of motion is the normal mode method. 1 In the normal 
mode method, the vector of a response quantity expressed in one coordinate 
system, such as the relative displacement between each node and the base { Z> , 
is transformed into the "normal" coordinate system by the relation 

{Z} = [<J>]{Y} (1.6) 

where 

[<J>] = a matrix of the natural mode shapes of the system with the 

ith column in the matrix representing the deflected shape 
of the ith mode of the system. 

and {Y} = a vector of modal response quantities with the ith element 

in the vector representing the modal response quantity in 
the ith mode. 

Thus, Eq. 1.6 expresses the response at each node in the model as a combina- 
tion of the responses at that node in each of the natural modes of vibration 
in the system. The natural mode shapes in the modal matrix [<}>] are orthogonal 1 , 
with respect to ["MJ and [K] and will be orthogonal with [C s ] if the elements 
°f [C ] are specified in a certain way. 21 Substitution of Eq. 1.6 into Eq. 1.5 
yields. 



6 



C'MJ[4>]{Y} + [C S ][*]{Y> + rK DL J[c(.]{Y} + [K][<j>]{Y} 

= -["Mj{x g } - rK DL J{x g } (1.7) 

Pre-multi plying Eq. 1.7 by yields, 

[*] T rMj[*]{Y} + [4 »] T [C s ][ 0,]{Y} + [<(»] T rK DL J[4»]{Y} 

+ [<J)] T [K][c()]{Y} = -U]YMJ{x g } - W T rK DL J(x g } (1.8) 

Recognizing the orthogonal relationships between the modal matrix [$] and 
other matrices in Eq. 1.8, the following matrices are defined: 



t~M*J = M T ['MJM 
rc*j = u] t [c s ][4»] 

[■K*J - M t [k]M 
[Kq L ] = C*] t C-k dl J[*] 

r-</j = 



a diagonal modal mass matrix 

a diagonal modal structural damping 
matrix. The form of [C s ] which permits 
this is discussed later, 
a diagonal modal stiffness matrix 

a modal linear hydraulic damping matrix 
which is not necessarily diagonal 
a matrix of circular natural frequencies 
of the system 



Using these quantites, Eq. 1.8 becomes 

N\j{Y} + rc*J{Y} + [Kg L ]{Y} + M<*-J{Y} 
= -ufr-Mjrxg} - [ ( f)] T r'K DL j{x g } 



(1.9) 
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If the modal linear hydraulic damping matrix [K^] can be made dia- 
gonal, then the set of simultaneous differential equations of motion repre- 
sented by Eq. 1.9 will be completely uncoupled. Suppose that the nondiagonal 

matrix [K^] can, in some satisfactory manner, be converted into a diagonal 
* 

matrix, say -]• Then Eq. 1.9 can be written 

[~M*J{Y} + + r K*^]{Y} 



•M T rMj{x g > - M T n< DL j<; g } 



( 1 . 10 ) 



Or, for the ith mode. 



•• 'k • 'k 

M.. Y. + C Y. + K n| Y. + K.. Y. 
n l s^. l DU.. i ni 



nn 




nn 


y <j).. m. . 

j=i J1 

/ 


(x ) - 
g 


l <t> • • K m 
0=1 DL oo | 



(x ) 

g 



where 

nn = number of nodes in the model 
Equation 1.11 can also be written as 



(i.n) 



Y i + 2s Sl “ii Y i + 2s h 1 “ii Y i + “ii Y i 



nn \ 

. 1=1 ** ^ 



★ 

M 11 



<V - 



I nn 

l 

jfl 



,L *ji K DL.. 



* 

M n 



ii 



/ 



<V 



( 1 . 12 ) 



where 



s ii 



s . 



1 2M ii“ii 



structural damping ratio in the ith mode 
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* 

K m. 

8. = — * = linear hydraulic damping ratio in the ith mode 

If the modal response quantity of interest, Y, Y, or Y can be ob- 
tained from Eq. 1.12 for each mode, then the total response of the structure 
can be obtained by using the relationship given in Eq. 1.6. The modal re- 
sponse quantities vary with time, and the time variation is different for 
each mode, so an exact calculation of the total response of the structure 
can only be made if the exact time histories of the modal responses are known 
and superimposed on one another. Generally, the maximum response quantity 
at each node of the model is desired, but the modal maxima do not necessarily 
occur at the same time. An approximation to the actual maximum response at 
each node of the model can be made by determining the maximum response in each 
mode and then combining the modal responses in some fashion. The response 
spectrum is a convenient way to determine the maximum modal response. 11-13 

1.3 Summary of Previous Research 

Equation 1.3, or a similar form, has been solved by other investi- 
gators using various methods. Much of the existing work in this area has 
been devoted to the response of offshore structures to wave forces, but the 
nonlinear difficulties are somewhat similar for wave forces and earthquake 
forces acting on submerged or partially submerged structures. Further, the 
bulk of recent investigations has been devoted to the use of stochastic repre- 
sentations of wave and earthquake motions. 

Burke and Tighe 4 integrate the equations of motion directly for a 
series of towers subjected to different deterministic ground acceleration 
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records by using an iterative technique which computes the response quanti- 
ties of interest at step-by-step time points throughout the earthquake record. 

For a tower in 800 feet of water subjected to four accelerograms (normalized 
to the intensity of the 1940 El Centro earthquake), the effect of the water- 
structure interaction was found to reduce the maximum base shear and the maxi- 
mum overturning moment by about 33 percent from the values found when the 
water-structure interaction was not included in the solution. This method, 
while permitting the direct use of the structural properties of the system, 
including the nonlinear drag force, requires the use of a digital computer 
and the availability of actual or simulated earthquake records. For purposes 
of preliminary design of a structure for a "design" earthquake of the designer's 
choice, however, the step-by-step technique is costly, time consuming, and 
requires that a time history of the "design" earthquake be available. Pre- 
liminary design would best be served by a rapid, approximate method of solu- 
tion which permits the use of selected parameters for the design earthquake 
and accounts for the effects of the nonlinear drag forces, yet yields satis- 
factory results. 

Penzien, in association with other investigators, 6 " 8 linearizes 
the drag forces in Eq. 1.3 by replacing the nonlinear force vector with a 
linear force vector which performs the same amount of work during the history 
of the response. This linearization results in an "equivalent" linear hydraulic 
damping matrix acting on the system in much the same manner as the structural 
damping matrix acts. The linear equations of motion are then transformed into 
a normal coordinate system and solved by the normal mode method of superposition. 1 
Unfortunately, the linear hydraulic damping matrix does not uncouple when the 






I 
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equations are transformed into the normal coordinate system. As a result, 
an uncoupling (or diagonalizing) technique, similar to the original linear- 
izing technique, must be employed in order to separate the equations into the 
desired modal equations. The linearizing and diagonalizing are accomplished 
by an iterative process which uses actual response quantities of the system 
as they are developed, and the procedure is time consuming and requires great 
care in its formulation and use. For towers subjected to a stochastic repre- 
sentation of the ground motion, equivalent linear hydraulic damping ratios on 
the order of 9 percent of critical are reported for the first two modes of 
towers in deep water, with "very much smaller" hydraulic damping ratios for 
the tower in shallow water. 6 

Foster 9 solves the problem of an offshore tower subjected to 
stochastic representation of wave forces by linearizing the drag force, then, 
recognizing that nonproportional hydraulic damping may result, solving the 
linear equations of motion using complex eigenvalues and eigenvectors. This 
procedure recognizes that, unless the linearization of the drag force is done 
in a certain way, classical normal modes may not necessarily exist in the equa- 
tions, thus the normal mode method of superposition is not a rigorously cor- 
rect approach. 

Nath and Harleman 10 conclude that, for large diameter cylinders in 
deep water waves, the nonlinear drag force can be neglected, and, since shallow 
water structures are relatively stiff, and therefore the dynamic response to 
waves is not greatly magnified over the static response, design of shallow water 
structures can proceed on a static basis. These conclusions, however, were 
arrived at by considering the response of structures to waves, which generally 
have much longer period components of motion than do earthquakes. 
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1.4 Objective and Scope 

The objective of this thesis is to investigate the effects of the 
nonlinear hydrodynamic drag force on the earthquake response of single degree 
of freedom and multi -degree of freedom systems and to develop a satisfactory 
method of linearizing this force in order that conventional methods of solving 
linear differential equations may be employed in the solution of the equations 
of motion. 

A further objective is to develop an approximate method for the 
solution of Eq. 1.3 by the superposition of normal mode responses, using 
response spectra to estimate the modal response quantities necessary. A rapid, 
yet fairly accurate method of solution is desirable in the preliminary stages 
of the design of an offshore structure in order that the designer may have 
an estimate of the response of a tentative design to one or more design earth- 
quakes of his specification. It should be kept in mind that the approximate 
solution developed in this thesis is no substitute for a more detailed analysis 
of the structure once the design is in its final stages. 

In Chapter 2, the effects of the nonlinear hydrodynamic drag force 
on the response of a single degree of freedom oscillator are explored. A 
reasonable method of linearizing the drag force and a method for conservatively 
predicting the resultant linear drag effects are developed and shown. A method 
is developed for estimating certain response quantities of the oscillator in 
order to determine more accurately the range of hydraulic drag and inertia 
coefficients which should be used in the force equation. The response of the 
oscillator in the inelastic region, assuming elasto-plastic force-displacement 
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characteristics of the spring, is shown. Finally, the construction of ideal- 
ized response spectra, which include the effects of hydraulic drag, is shown. 

The concepts developed in Chapter 2 for the single degree of freedom 
system are extended to the multi -degree of freedom system in Chapter 3. The 
nonlinear drag force is linearized in such a manner as to insure the existence 
of classical normal modes in the equations of motion in order that normal mode 
superposition may be used as an approximate solution to the problem. The method 
of linearization is shown to be accurate for a series of five degree of free- 
dom systems which are representative of the type of structure being investi- 
gated. A means of estimating the amount of equivalent linear hydraulic damping, 
using ground motion velocity characteristics, is developed and is shown to 
yield good results. The construction of an idealized response spectrum for 
each mode of the multi -degree of freedom system is discussed and described, 
and sample problems are solved by the superposition of normal mode responses. 

The application of the approximate method to a series of actual off- 
shore towers is demonstrated, and the results of calculations of the response 
of these towers are discussed in Chapter 4. 

1.5 Related Topics 

In this section, several topics of general interest in the field of 
earthquake engineering are briefly discussed in order that the reader may more 
fully understand the display of results by the use of response spectra, the 
numerical method of integrating the equations of motion, and the characteris- 
tics of the earthquake records used in the research. 
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1.5.1 Response Spectra 

The equation of motion for a damped, single degree of freedom os- 
cillator in air, subjected to a base disturbance is 

mx s + C s (x s - x ) + K(x s - x g ) = 0 (1.13) 

By making the transformation Z = x $ - x g , Eq. 1.13 becomes 

mZ + C Z + KZ = - mx (1.14) 

^ y 

Or, in an alternate form 

Z + 2e uZ + u 2 Z = - x (1.15) 

s 9 

The solution of Eq. 1.15 depends on the characteristics of the ground motion 
function, the amount of damping in the system as represented by the dimension- 
less parameter B s > and the natural frequency of the system, w. The frequency 
parameter enters into the solution as a dimensionless product of w and a 
characteristic time of the base motion. 13 

A plot of the maximum value of a response quantity, without regard 
to sign, as computed from Eqs. 1.13, 1.14 or 1.15, as a function of the fre- 
quency of the system is called a response spectrum for that quantity. The 
response spectrum is a convenient way to represent various maximum response 
quantities of interest as a function of the variables involved in the solution. 
The response quantity of greatest interest is usually the maximum spring dis- 
tortion, |Zl . which occurs during the response history of the oscillator, 
max 

although other response quantities are sometimes of interest. 13 

The response spectrum is conveniently shown on tripartite logarithmic 
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graph paper as shown in Fig. 1.2. The grid of this representation is con- 
structed such that the abscissa represents the frequency of the oscillator, 
the ordinate represents velocity units of response, the diagonal lines 
sloping up to the right are lines of constant displacement units, and the 
lines sloping down to the right are lines of constant acceleration units. 

In this thesis, unless otherwise indicated, the response spectra shown are 
deformation response spectra in which the pseudo-velocity, oj | Z | is plotted 
against the undamped frequency of the system. The displacement quantity 
shown is the maximum spring deformation, I Z I , and the acceleration quantity 

ITIaX 

o 

is the pseudo-acceleration, to |Z| . Response spectra for a damped single 

degree of freedom oscillator, in air, are shown in Fig. 1.2 for the El Centro 
earthquake. The shape of the spectra is typical for strong motion earthquakes 
of the type used in this thesis. More detailed information on response spectra 
can be found in the literature. 11-14 

Since the response spectra for strong motion earthquakes exhibit the 
characteristic shape shown in Fig. 1.2, the spectra can be idealized by straight 
line segments. Amplification factors, which relate the idealized spectrum to 
the ground motion quantities x^ max » x^ max , and x^ max , have been developed 
for a range of damping ratios. 14 The ability to construct an idealized repre- 
sentation of the response spectrum for a given set of conditions permits an 
estimate of the response to be made without having an actual earthquake record 
available. 

1.5.2 Numerical Integration of the Equations of Motion 



The equations of motion for a system subjected to a base disturbance 
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can be solved by numerical integration. In this thesis, numerical integra- 
tion of the equations of motion is performed by the Newmark "Beta Method." 5 
In the calculations, $ is set to be equal to 1/6, thus representing the ground 
response acceleration as linear between adjacent time points. The time inter- 
val used is 1/20 of the natural period for frequencies up to 10 cps and 1/40 
of the natural period for frequencies equal to, or greater than 10 cps. These 
time increments insure good accuracy, stability and convergence in the solu-' 
tion. Convergence of the solution at any time point is assumed to have oc- 
curred if two successive iterations yield values of acceleration which are 
within 0.1 percent. In the case of multi -degree of freedom equations of 
motion, the time increment is based on the shortest natural period in the 
system, and iteration is carried out until the computed acceleration at each 
node is within 0.1 percent of the value computed on the previous cycle. In 
all cases, the integration of the equations of motion is carried out for a 
period of time equal to one-half the natural period (the longest natural 
period in the case of multi-degree of freedom systems) after the base distur- „ 
bance ends. 

1.5.3 Earthquake Records Considered 

Four actual earthquake records are used in this thesis. The records 
are described numerically in Table 1.1 and the time history of the ground ac- 
celeration, velocity, and displacement for each record is shown in Figs. 1.3 
through 1.6. In each earthquake record, a very small (less than .005g) arti- 
ficial acceleration point has been generated near the end of the record so 
that the velocity of the base motion is zero at the end of the record. 
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Although a parabolic baseline correction 16 had previously been applied to the 
acceleration records, a slight residual velocity still remained at the end of 
the record and it is this residual velocity which is removed by the artificial 
acceleration point. The four records chosen are considered to be representa- 
tive of strong motion earthquakes on firm soil. 

1.6 Notation 



degree of freedom oscillator and the multi-degree of freedom system, with 
matrix and vector notation indicating the multi -degree of freedom system. 
Symbols are defined as they first appear in the text. For ease of reference 
the symbols used are summarized in this section. 



As much as possible, the same symbols are used for both the single 



r-Aj 



projected area matrix 



C 



D 



empirical dimensionless drag coefficient 



C 



I 



empirical dimensionless inertia coefficient 



[C s ] 



structural damping coefficient matrix 




modal damping coefficient matrix = [<f>]^"[C s ][<j>] 



D 



diameter of a single degree of freedom oscillator 



{HDPF} 



vector of participation factors for hydrodynamic 
parameter response quantities 



{HDPF} = ^K*^J _1 [4.] T [K]{I} 



{ 1 } 



unit vector 
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[K] 


structure stiffness matrix 




modal stiffness matrix = [4>]^[K][<J>] 


r-K D j 


drag factor matrix = ^-pC'CpAJ 


«£> 


modal drag factor vector for very high frequency 
system = ["M [<1>] T E' KpJ{I} 


rK DL ~J 


linear drag factor matrix 


i — i 

_j 

<? Q 

i i 


modal linear drag factor matrix = [<j>]^[~ Kp L ~J[<j>] 


K n x 

D q max 
K 


dimensionless drag parameter 


t'M-J 


1 

total structure dynamic mass matrix = ['mJ + ['m^J 




modal mass matrix = [<|)]^["MO[(J)] 


{p> 


hydrodynamic force vector 


{PF} 


participation factor vector for displacement response 

= r'M 




enclosed volume matrix 


{Y}, {Y}, {Y} 
{Q}, {Q>, {Q> 


displacement, velocity , and acceleration vectors 
in normal coordinates 


(Z>, {Z}, {Z} 


vectors of displacement, velocity, and acceleration 
of the structure relative to the base 


Z 

y 


deformation at which yielding occurs 
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a 


relative acceleration between fluid and a single 
degree of freedom oscillator 


i 


subscripts indicating mode number 


j 


subscripts indicating node number 




structural mass matrix 


[Vj 


added mass matrix = pC'CjV.J 


nm 


number of modes in the multi -degree of freedom system 


nn 


number of nodes in the lumped parameter model 


P(t) 


hydrodynamic force acting on the single degree of 
freedom oscillator 


V 


relative velocity between fluid and a single degree 
of freedom oscillator 


{X gl, (X g }, (X g ) 


ground displacement, velocity and acceleration vectors 


• 

x g max’ x g max’ x g max 


maximum displacement, velocity, and acceleration in 
the ground motion record 


{x s h {i s }, {x s } 


vectors of structure displacement, velocity, and 
acceleration 


y, , y, , y,. 
V V “i 


displacement, velocity, and acceleration of the ith 
mode single degree of freedom equation of motion 
after the modal participation factor has been taken 
out. Expressed in normal coordinates 
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a 




{y} 

V 



[♦] 

P 



“in- 



constant of proportionality between the linear drag 
factor matrix and the total mass matrix 
= a[~ M 



linear hydraulic damping ratio in the ith mode 

structural damping ratio in the ith mode 
C 

6 S = V 1 — 

vector of modal stiffness factors {y} = [<f>]^[K]{I) 



ductility factor for inelastic response 

IZI, 



'max 



mode shape matrix 



mass density of the fluid 



circular natural frequency of the ith mode 

l> 2 j = [Vj'Vif-J 
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CHAPTER 2 

SINGLE DEGREE OF FREEDOM OSCILLATOR 



2.1 General 

The normal mode method of solving Eq. 1.3 requires the solution 
of single degree of freedom equations of motion. For this reason, and to 
understand more fully the effects of the nonlinear hydrodynamic drag force 
on the response of a simple system. Chapter 2 is devoted to the investigation 
of the effects of the hydrodynamic drag force on the response of a single de- 
gree of freedom oscillator. 

In the investigation of the response of the oscillator with non- 
linear drag characteristics it is shown that the response of the oscillator 
in the very high frequency region can be estimated if the drag properties of 
the system are known and if the maximum velocity of the base motion is known. 
The nonlinear drag force is linearized by equating the work done by a ficti- 
tious linear force with the work done by the real, nonlinear force. This 
method of linearization is shown to be accurate until the high frequency re- 
gion of the spectrum is reached. It is shown that the linear drag force is 
equivalent to a linear hydraulic damper in which the force is dependent upon 
the velocity of the mass, and that the magnitude of the hydraulic damping 
coefficient is dependent upon the amount of structural damping in the system. 
A means of estimating the amount of linear hydraulic damping is developed 
for earthquake motions, and the estimate is shown to be conservative except 
in the low frequency range. The maximum response of the system in the mid- 
frequency range to earthquake excitation is shown to be dependent upon the 
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total amount of linear damping present, and not upon the ratio in which the 
damping is distributed between structural and hydraulic. 



Response quantities which may be used in order to select reasonable 



values of the inertia and drag coefficients, Cj and Cp, are explored, and a 
method of estimating these response quantities is developed. The inelastic 
response of an oscillator with elasto-plastic force-displacement character- 
istics of the spring is investigated and the results are shown in the form of 
response spectra. 



The construction of idealized elastic response spectra which include 



the effects of the nonlinear hydrodynamic forces is demonstrated in this 
chapter. 

2.2 Elastic Deformation Response Spectra 

The equation of motion can be written for the single degree of 
freedom oscillator as 



MZ + C Z + KZ 
s 




( 2 . 1 ) 



or 





K, 



( 2 . 2 ) 



or 





(2.3) 



As pointed out in Chapter 1, the solution of the equation of motion for a 
land-based single degree of freedom oscillator (Eq. 2.1 without the nonlinear 
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hydrodynamic drag force k q^ s I I ) depends on the characteristics of the ground 
motion, the amount of damping in the system as represented by the dimensionless 
parameter $ s , and the natural frequency of the system, u>. The solution of 
Eq. 2.1 depends on the magnitude of Kp in addition to these parameters. A 
dimensional analysis of Eq. 2.1 yields the dimensionless drag parameter 

K P *q max 
K 

as one of the parameters on which the response of the oscillator depends. 

When using a maximum ground acceleration of 0.5 g, values of 0.01, 0.2, and 

0.6 represent a reasonable range over which to vary the dimensionless drag 

parameter. For the maximum ground acceleration in the El Centro earthquake 

record (0.312 g), values of 0.006, 0.08, and 0.50 for the parameter yield 

slightly greater variations in the magnitude of Kp. Obviously, a zero value 

of the drag parameter represents a "no-drag" system and increasing values of 

the parameter represent increasing magnitudes of Kp in the system for fixed 

values of K and x . 

g max 

Equation 2.1 is solved numerically for the steady-state maximum 
response of the oscillator when the base acceleration is 

Xg = 16.1 sin 2irt 

where 

Xg is in ft/sec 2 
t is in seconds 

Deformation response spectra (in which z I max is plotted against frequency) 
are plotted in Figs. 2.1 and 2.2 to show the results of the calculations. 
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In order to investigate the response of the oscillator to a slightly 
more complicated representation of the base disturbance, the base acceleration 
is represented as a six term cosine series. The numerical characteristics of 
the series used in the response calculations are shown in Table 2.1. The 
series has the features of zero acceleration, velocity, and displacement at 
the beginning and end of the base motion. Equation 2.1 is solved numerically 
for the maximum response of the oscillator when subjected to five cycles of 
this base disturbance, and the results are shown in the form of deformation 
response spectra in Figs. 2.3 and 2.4. 

Equation 2.1 is solved numerically for the maximum response of the 
oscillator when subjected to the El Centro earthquake record. The results 
are shown in the form of deformation response spectra in Figs. 2.5 through 
2.7. 

A study of the spectra presented in Figs. 2.1 through 2.7 reveals 
several interesting characteristics which are best discussed in terms of 
frequency ranges of the oscillator. 

In the low frequency region, the maximum spring deformation, |Z| , 

approaches that for the case of K^x^ max /K = 0, regardless of the amount of 
hydrodynamic drag in the system for the range of drag parameter used. This 
is reasonable to expect, since the magnitude of the drag force is a function 
of the oscillator velocity, x $ , and very soft (low frequency) systems exhibit 
little velocity response of the mass (x s ). 

In the mid-frequency region, where spectra for the conventional 
"no-drag" oscillator show an amplification of the pseudo-velocity over the 
maximum ground velocity, the spectra for systems with drag also show an ampli- 
fication of the pseudo-velocity over the maximum ground velocity. However, 
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the presence of the nonlinear hydrodynamic drag force causes a reduction in 
the response of the oscillator from the response of the "no-drag" oscillator. 
The amount of reduction is greater for greater values of the drag parameter 
up to some frequency. This reduction would suggest that, in the mid-frequency 
range of the spectrum, the nonlinear hydrodynamic force can be replaced by 
an equivalent linear hydrodynamic force which will effect the response in 
much the same manner as the presence of structural damping effects the re- 
sponse. The amount of reduction in response caused by the drag force is also 
seen from the spectra to be dependent upon the amount of structural damping 
in the system. This is easily noted by comparing Figs. 2.5 and 2.6. In Fig. 
2.5, the system has no structural damping, and a slight amount of drag (repre- 
sented by KpXg max /K = 0.006) is seen to reduce significantly the maximum 
response of the oscillator in the frequency range of about 1 to 10 cps. In 
Fig 2.6 the system has 2 percent structural damping, and the drag parameter 
of 0.006 is seen to cause very little reduction in the response spectrum be- 
tween 1 and 10 cps. The effects of structural damping on the amount of re- 
duction in response caused by the hydrodynamic drag are to be expected, since 
the drag force is a function of the oscillator velocity, x s , and this velocity 
is a function of the amount of structural damping in the system. 

In the high frequency region, response spectra for "no-drag" oscil- 
lators typically show the pseudo-acceleration approaching the maximum ground 
acceleration, x^ max> Response spectra in the high frequency region for the 
oscillators with nonlinear hydrodynamic force exhibit entirely different 
characteristics, however. In all the response spectra shown, a frequency is 
reached at which the maximum response quantity |Z| max for the oscillator with 
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drag exceeds that for the "no-drag" oscillator. In fact, for the systems 
with larger values of the drag parameter, the quantity |Z| is seen to ap- 

ilia A 

proach a relatively constant value. The magnitude of this high frequency 
"tail" in the spectrum can be estimated as follows: 

Consider an undamped system with spring stiffness so great that 



X„ * x „ 

s g 



x„ * x 

s g 



The forces on the oscillator are 



F = F. .. + F, 

spring inertia drag 



^spring = M s + VA 1 



or. 



F «Mx ± K ( x Y 

spring g V X g' 



Since the drag parameter 



K n x 

D q max 
K 



is held constant, must increase as the spring stiffness K increases, so 
a frequency will be reached at which the inertia force can be neglected with 
respect to the drag force. If the inertia force is neglected, then 



^spring^max 



« Mx a y 

D g max 



And 



'max 



(F ) K 
v spring y max D 

K K 



(i g m ax>‘ 



1 / •• c 

K n x x„ , 

^ D q max ^ q max 



K 



) (2.4) 



g max 
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A comparison of |Z| obtained at high frequencies by the numerical 

max 

integration of Eq. 2.1 with |Z| computed from Eq. 2.4 is shown in Table 

I II Cl A 

2.2 for all systems and base disturbances considered. A study of the data 
presented in Table 2.2 shows that, particularly for the earthquake records, 
the high frequency tail of the deformation response spectrum can be satisfac- 
torily estimated by the use of Eq. 2.4. 

When an acceleration record is integrated numerically in order to 
obtain the velocity and displacement records, the integration is done from 
time point to time point in the record. This double integration yields the 

values of x = 12.79 in. /sec and x „ = 8.24 inches for the El Centro 

g max g max 

record as shown in Table 1.1. However, if the acceleration record is inte- 
grated to find the ground velocity and displacement at some time point which 
is between time points in the acceleration record, linear interpolation of 
the acceleration record must be performed in order to obtain the acceleration 
quantity necessary to compute the ground velocity and displacement at the 
intermediate time point. This interpolation introduces small errors into 
the acceleration, which in turn introduces slightly larger errors into the 
ground velocity which in turn introduces even larger errors into the ground 
displacement computed at the intermediate time point. 16 It was found that, 
when using the very small time increments necessary for stability and con- 
vergence in the numerical integration of the equation of motion at high fre- 
quencies (above about 8 cps), the interpolation of the El Centro acceleration 

record introduced errors into the magnitude of x such that x was found 

g g max 

to be about 14.0 in. /sec in lieu of the value 12.79 in. /sec previously ob- 
tained, even though double precision was used in all calculations. In order 
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to achieve comparability of results in the comparison of |Z| estimated from 

1 1 IQ A 

Eq. 2.4 and I z I max obtained at high frequencies from the integration of the 



solve Eq. 1.3 depends on the ability to linearize the nonlinear hydrodynamic 
force vector Kp { x s | x s | } in a satisfactory manner. The linearization of 
the hydrodynamic drag force in the single degree of freedom system is inves- 
tigated in this section. It is desirable to replace the nonlinear drag 
force, Kp^ s l^ s l» ln Eq. 2.1 by a linear drag force, say, K^x , which will 
yield the same maximum response in the system. Consider the work done by the 
two forces: 15 



equation of motion, x =14.0 in. /sec is used in Eq. 2.4 to obtain the 

CJ ITIaX 

results shown in Table 2.2 for K n 5i _ /K = 0.006 and 0.08. 

u 9 max 



2.3 Linearization of the Drag Force 



It was seen in Chapter 1 that the use of the normal mode method to 



Work done by the nonlinear 
drag force over a period of time 



(V s |x s l )( x s )dt 



0 



t 



Work done by a linear drag 
force over a period of time 



(K DL ; s )( x s ) dt 



0 



Equating the amounts of work done, 



o 




o 






(2.5) 
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(The symbol < > denotes the average value of the quantity enclosed over the 
time of the response of the oscillator.) This expression is the same expres- 
sion used by other investigators 6-9 to linearize a nonlinear force. With the 
nonlinear drag force replaced by an "equivalent" linear drag force, Eq. 2.1 
becomes , 

MZ + C Z + KZ = -Mx - K nl x 

s g DL s 

MZ + C g Z + K dl Z + KZ = -Mx g - K DL x g (2.6) 

or, 

Z + 26 s wZ + 26 h wZ + o) 2 Z = -x - 26 h wx g (2.7) 

where, 

^DL 

^h = 2Mw = ec l u ‘' va 1 ent linear hydraulic damping ratio 

It now remains to investigate how the maximum response in the equi- 
valent linear system compares with the maximum response in the actual nonlinear 
system. 

Equation 2.6 is solved numerically for |Z| max using a linear drag 
factor Kp^ which is linearized in accordance with Eq. 2.5 for the single de- 
gree of freedom oscillator subjected to the series representation of the base 
disturbance. The quantities < x^|x s |> and <x^> necessary for the linearization 
were obtained from the integration of the nonlinear equation of motion for the 
same set of conditions. The results are plotted on Figs. 2.8 and 2.9, and it 
is seen that, for this base disturbance, the linearization is very good until 
the high frequency tail begins to appear. This might have been expected. 
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since a linearization of a weak nonlinearity would normally be better than a 
linearization of a strong nonlinearity, and the high frequency tail is the 
result of the nonlinear hydrodynamic drag force becoming quite strong and 
dominant in Eq. 2.1. 

A similar set of calculations is made for the El Centro earthquake. 
The response spectra for Kp and Kp^ are shown in Fig. 2.10 for = 0.02 and 
KpXg max /K = 0.50. It is seen from this figure that the linearization for 
the earthquake motion, while not quite as good as the linearization for the 
series representation of the base motion, is still reasonably good up until 
the high frequency tail of the spectrum is reached. A comparison of |Z| max 
computed for the nonlinear drag system with |Z| for the linear drag sys- 
tern is shown in Figs. 2.11 through 2.13 for all of the variables considered. 
It is seen from a study of these figures that the maximum response of the 
linear system is close to the maximum response of the nonlinear system until 
the high frequency region is reached. The figures further indicate that the 
linearized response is slightly on the conservative side of the actual re- 
sponse until the high frequency region is reached. The accuracy of the 
linearization in the high frequency region is expected to become poor as the 
nonlinearity is becoming more and more dominant in the system, and thus the 
linearization becomes less and less equivalent. However, Eq. 2.4 permits a 
fairly good estimate to be made of the maximum relative displacement in the 
high frequency region. 

Since the linearization of the drag force in Eq. 2.1 is seen to 
yield a slightly conservative maximum response at frequencies below the be- 
ginning of the high frequency tail in the response spectrum, a method of 
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estimating the linear drag factor from the known nonlinear drag factor 
Kp is now developed for earthquake motions. A study of Eq. 2.5 shows that 
two velocity quantities, <x^ | x s j > and <x4 are necessary in order to linearize 
the drag force and, since these quantities are actual response quantities 
of the nonlinear system, they cannot be calculated a priori . A series of 
computations is macfe in which the quantity 

< >'sl x sl > 

4 

is computed for the actual, nonlinear oscillator subjected to the El Centro 
earthquake. The results are shown in Figs. 2.14 through 2.16. From these 
figures it is seen that, for the range of variables studied, if the quantity 



• L i • i 

<X S I* S I> 

4 



is estimated as 

•2 1 • . 

<x x > 

4 

for the particular earthquake considered, then this estimated value will be 
less than the actual 



•ti • I 

<x„ x > 
s' s' 

4 



for all frequencies except low frequencies. Using the value of 
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in Eq. 2.5 to linearize the drag term results in a value of which is less 
than the value of which would result from using the actual oscillator 
velocity quantities. When this lower value of is converted into the 
linear hydraulic damping ratio, shown in Eq. 2.7, the value of 3^ is 
less than the actual value of 3 h which would be computed had the actual os- 
cillator velocity terms been available for use in Eq. 2.5. In other words, 
linearizing the drag term by using 



•2 1 • i 

< Xql X ql > 

•2 

<x g > 



in Eq. 2.5 leads to a conservative value of which, in turn, leads to a 
conservative value of the maximum response in all except the low frequency 
region and the high frequency tail region of the spectrum. The unconserva- 
tive nature of K QL in the low frequency region when computed from the ground 
velocity terms is reasonable to expect, however, as the amount of hydraulic 
damping is dependent on the velocity response of the oscillator, and very 
soft (low frequency) systems have little velocity response, so the amount of 
hydraulic damping should approach zero as the system becomes softer and softer. 
Figures 2.14 through 2.16 further show that, for a system with given drag 
characteristics, the linear drag factor estimated by using the quantity 



• 2 1 • i 

<xlx > 

g'V 

4 
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becomes less conservative as the amount of structural damping increases and 
as the amount of nonlinear drag increases. The linear drag factor estimated 
in this manner does not become unconservative (except in the low frequency 
region) for the range of variables studied. 

The amount of conservatism introduced by the use of 



<x q ' x a' > 

4 

in linearizing the drag factor, Kp, is seen in Figs. 2.17 through 2.19. In 
these figures, the spectrum for the linear drag solution is seen to be con- 
servative up until the high frequency tail begins. The linearized spectra 
were found to be slightly unconservative, as expected, in the low frequency 
region, but the scale of the figures is such that this unconservatism cannot 
easily be seen. The use of Eq. 2.4, repeated here for convenience. 



l z U ■ <T><W 2 ■ 

3 g max 



(2.4) 



is also seen to provide a very good estimate of the nonlinear drag spectrum 
in the high frequency region. At the point in the spectrum at which the high 
frequency tail just begins to appear, the intersection of the linear drag 
spectrum with the estimate of the high frequency tail provided by Eq. 2.4 
is seen to be slightly unconservative. The amount of unconservatism in this 
transition "saddle" to the high frequency tail is reduced as the amount of 
nonlinear drag in the system increases. It is seen from the figures that, 
while linearizing the drag force by the relation, 

<x _ I x | > 

k = k ( 9 _g \ 

K 0I_ 1 



4 



( 2 . 8 ) 
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yields conservative results in the mid-frequency region of the spectrum, the 
degree of conservatism is not so great as to be considered too conservative. 

The presence of the nonlinear hydrodynamic drag force has been seen 

to cause a reduction in I Z I from that computed for the "no-drag" oscillator 

' 'max 

at a given frequency, up until the high frequency tail region is reached. A 
measure of how much of this reduction is accounted for in the computation of 
| Z | for the system with linear drag (linearized by the ground velocity) 

IliuX 

is seen in Fig. 2.20. The ordinate in the figure is the ratio of the reduc- 
tion (from the "no-drag" case) in I Z | max when the actual nonlinear drag force 
is used to the reduction (again, from the "no-drag" case) in I z I max when tbe 
drag force has been linearized by the ground velocity terms. In Fig. 2.20, 
values of the ordinate greater than 1 represent unconservative results. The 
unconservative nature of the linearization in the low frequency region and 
the high frequency tail region is noted. While the linearization of low-drag 
systems is seen to account for less than 50 percent of the actual reduction 
in l z L;w in man y cases, the actual reduction for low-drag systems is quite 
small itself, as seen in Fig. 2.17. It is further seen from the figure that 
more of the reduction in |Z| max is accounted for in the high-drag systems than 
in the low-drag systems. 

Figures 2.17 through 2.20 show a measure of the amount of conserva- 
tism introduced in I Z | max by using the ground velocity characteristics to 
linearize the drag force when the structural damping ratio, 3 s , is 0.02. 

For systems with less structural damping, the amount of conservatism will 
be slightly greater, and for systems with more structural damping the amount 
of conservatism will be slightly less. 
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2.4 Linear Hydraulic Damping 

The preceding section of this chapter has demonstrated that the 
drag factor can be linearized in a manner which yields slightly conserva- 
tive values of |Z| max over a wide range of frequencies in the response 
spectrum. The linearized drag factor, K^, can be converted into a linear 
hydraulic damping ratio, 3^> as shown by Eq. 2.7, which is repeated here for 
convenience. 

Z + 2B s toZ + 23 h o)Z + oo 2 Z = -x - 2B h wx g (2.7) 

Equation 2.7 is seen to have the same form as Eq. 1.12 with the exception of 
the coefficients of the forces on the right side of the equation. Since the 
solution of Eq. 1.12 was shown in Chapter 1 to be necessary in order to use 
the normal mode method of solution for the response of the multi -degree of 
freedom model, the effects of B h on the solution of Eq. 2.7 will now be in- 
vestigated. 

Equation 2.7 is solved numerically for |Z| „ for all four earth- 

quakes with values of B $ ranging from 0 to 0.10 and values of 3^ ranging from 
0 to 0.10. These ranges of B s and B^ are considered to be inclusive of all 
reasonable values of these variables for the type of structure considered in 
this thesis. The results of the computations are shown as response spectra 
in Figs. 2.21 through 2.24. These response spectra are seen to have the same 
general shape as spectra for an oscillator without hydraulic damping until 
the high frequency region is reached. It should be noted that a high fre- 
quency tail also appears in the spectra for systems with equivalent linear 
hydraulic damping. The magnitude of this tail can be shown to approach 
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20 . . 

( — -)(x ) from above. However, since the equivalent linearization of the 

nonlinear hydrodynamic drag force has already been shown to break down as 
the drag force becomes more dominant in the solution at high frequencies, an 
estimate of I Z | in the high frequency tail region of the spectrum for a 
given system can best be made by the use of Eq. 2.4. 

Newmark et al. 11,13 » 14 construct an idealized representation of 
the single degree of freedom response spectrum for an oscillator without drag 
forces by using amplication factors for displacement, velocity, and acceler- 
ation. These amplication factors are multiplied by the corresponding maximum 
ground motion quantity in order to obtain an idealized spectrum for an earth- 
quake motion of known (or assumed) values of x , x , and x . If 

g max g max g max 

amplification factors can be developed for the response spectrum of the system 
with equivalent linear hydraulic damping when 0 s and 0^ are known, and the 
high frequency tail of the spectrum for the actual nonlinear system is esti- 
mated from Eq. 2.4, then an idealized response spectrum can be constructed for 
a given set of conditions for the oscillator with nonlinear hydrodynamic drag 
forces. 



During the course of establishing the amplification factors necessary 
for the construction of idealized response spectra for systems with linear 
hydraulic damping ratios (Eq. 2.7), it was found that, within the range of 
variables studied, the amount of total linear damping, 0 $ + 0^, is the impor- 
tant factor controlling the response quantity |Z| throughout the spectrum 
except for the high frequency region. Comparisons of |Z| max computed from 
Eq. 2.7 at 31 frequencies from 0.04 cps to 15 cps were made for all four 
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earthquake records for total damping ratios equal to .02, .05, .07, .10, .12, 
and .15 with the damping being distributed between structural and hydraulic 
as shown in Table 2.3. The comparison showed that, up to the beginning of 
the high frequency tail in the spectrum, for a given value of the total damping 
ratio, IZl computed for the first distribution of 3 and g, was within 5 
percent of I Z I „„ computed for the second distribution of g and g, . Further, 
the great majority of the comparisons revealed a difference of less than 2 
percent between I ^ I max computed for the two different distributions of damping 
at a given value of the total damping ratio. The ratio 



m 

m 



max 



max 



for g $ = A, g h = B 

for g = B, g = A 
s h 



was found to be sometimes greater than 1 and sometimes less than 1, indicating 
that at least a small part of the difference may be the result of the conver- 
gence criterion (0.1 percent) used in the numerical integration procedure. 

A comparison of I Z I , w when g =0.00 and g, = 0.15 with |Z| , when 
g s = 0.15 and g^ = 0.00 was made for each of the four earthquakes considered. 
It was found that, for frequencies up to the beginning of the high frequency 
tail of the spectrum, the value of I Z | max computed for the first distribution 
of total damping was not always within 5 percent of the value of |Z| max com- 
puted for the second distribution of total damping. Most of the comparisons 
differed by less than 10 percent, but a few showed differences on the order 
of 15 percent. 

The conclusion to be drawn is that, for frequencies up to the be- 
ginning of the high frequency tail region, |Z| computed from Eq. 2.7 is 

I lid A 
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dependent upon the total linear damping in the system, and not the distribution 
of damping between structural damping and equivalent linear hydraulic damping 
for values of 3^ up to and including 0.10. Thus, the use of the amplification 
factors published in Ref. 14, and included in this thesis as Table 2.4, is a 
valid and proper method of constructing a portion of an idealized response 
spectrum for the single degree of freedom oscillator with drag characteristics. 

The scale of the response spectra shown in Figs. 2.21 through 2.24 
is so small that the conclusion is difficult to see from the figures. For 
illustrative purposes, a comparison of response spectra for 3 = .07 is shown 
in Fig. 2.25 for the Ferndale earthquake. It is readily seen from this figure 
that, up until the high frequency tail begins to appear, the spectral values 
are a function of the total amount of linear damping in the system, and the 
relative distribution of the total damping between structural and hydraulic 
has little effect. 

2.5 Idealized Elastic Deformation Response Spectrum 

It was mentioned in Section 1.5.1 that the characteristic shape of 
response spectra for earthquake motions has led to the development of amplifi- 
cation factors from which an idealized response spectrum can be constructed 
for a given set of conditions. The idealized spectrum permits an estimate 
to be made of the maximum response of the system. This response quantity 
may then be used as a modal response quantity for a particular mode of a 
multi -degree of freedom system, and the various modal response quantities 
may then be combined in some manner to yield an approximate solution for the 
multi-degree of freedom system as outlined in Chapter 1. 

The principles and concepts developed in Sections 2.2 through 2.4 
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are summarized here for the reader in such a manner as to demonstrate the 
construction of an idealized response spectrum for a single degree of freedom 
system, with nonlinear hydrodynamic drag forces, subjected to a strong-motion 
earthquake of the type used in this thesis. 

To construct an idealized deformation response spectrum for a given 



system: 



1. For the system, compute the drag factor, K Q , and the undamped 
circular natural frequency of the system, w. 

2. Establish the structural damping ratio, 8 , to be used in the 
system. For the type of structure discussed in this thesis, 

a structural damping ratio of between 0.005 (when the stress 
is below one-quarter of the yield stress) and 0.05 (when the 
stress is at, or just below, the yield stress) is reasonable 
to assume. 14 

3. Select the values of x„ m . x„ m . x„ m3x/ for which the 

g max g max g max 

spectrum is to be constructed. 

4. Compute the quantity 



K d ( 



• 2 1 • i 

<X X > 

£l g ) 

•2 ' 

<x > 

q 

2 M u 



It should be noted from Table 1.1 that the magnitude of 



<X ql X ql > 

•2 

<X 9 > 



is approximately one-half the magnitude of x for the four 

g max 

earthquakes shown. 
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5. Enter Table 2.4 with the value of expressed as a percent, 
obtain the displacement amplification factor, multiply the 
displacement amplification factor times max , and plot the 
result on tripartite logarithmic graph paper as a displacement 
line. This is shown as line 1 on Fig. 2.26. The method of 
linearizing the drag force by using ground velocity character- 
istics has been shown to be unconservative in the low frequency 
region. The use of 6 alone means that any hydraulic damping 
effects in the low frequency region are neglected. This is 
reasonable, however, because the small amount of velocity re- 
sponse in the low frequency region results in minor amounts of 
hydraulic damping. Further, as the frequency of the system 
becomes lower and lower, the amount of damping of any type has 
little effect on the maximum deformation. 

6. Enter Table 2.4 with the amount of total damping (3 g + as 
computed above) and obtain velocity and acceleration amplifi- 
cation factors. Multiply the velocity amplification factor 
times Xg max and plot the result as a horizontal line. Multiply 

the acceleration amplification factor times x and plot the 

g max r 

result as an acceleration line. These two lines are shown as 
lines 2 and 3 on Fig. 2.26. 

7. Plot the quantity (~j^)(Xg max J as a displacement line. This 
line is shown as line 4 on Fig. 2.26. 

8. The envelope formed by the intersection of the four lines plotted 
is the idealized response spectrum for the system and conditions 
given. 
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With the idealized response spectrum constructed as described in 
(1) through (8) above, an estimate of the maximum response of the system is 
obtained by entering the figure with the natural frequency of the system and 
reading |Z| m . pseudo-velocity (w|Z| ), or pseudo-acceleration (w |Zl ). 

If the characteristics of the system are changed, a new response spectrum 
should be constructed. 

Linear interpolation can be used to obtain amplification factors 
for values of damping not shown exactly in the table. 

2.6 Hydrodynamic Parameters 

The hydrodynamic force exerted by the fluid on the mass of the oscil- 
lator as it moves through the fluid is represented by Mori son's formula, Eq. 1.2 
as described in Section 1.2. This force expression contains two dimensionless, 
empirical coefficients, Cj and Cp, which are dependent upon the response of 
the oscillator for their magnitude. Many investigators have undertaken 
studies to determine the value of these coefficients under various conditions. 
The results of some of these investigations are published and discussed in 
Refs. 17 and 18. Morison's investigations 2 show a variation of the drag co- 
efficient Cp with Reynolds Number, which is related to the relative velocity 
between the mass and the fluid. Keulegan and Carpenter 19 have shown a variation 
in both coefficients with a dimensionless period parameter which can be in- 
terpreted as the ratio of the relative displacement between the mass and the 
fluid to the pile diameter for forces against circular piles. Laird 20 has 
shown that the two coefficients may vary greatly for flexible oscillating 
cylinders. Evaluation of the two coefficients for oscillatory motion is further 
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complicated by turbulence associated with eddy formations on either side of 
the mass as it moves back and forth through the fluid. 18 This thesis con- 

I 

siders that the drag factor, K^, and the added mass of the oscillator, m , 
remain constant throughout the entire time history of the oscillator motion. 

I 

The values of K Q and m used in the calculations in this thesis do not pre- 
suppose any specific values for Cj and C Q , except as values were assigned 
to these coefficients by Burke 4 in the formulation of properties of a series 
of offshore towers which are used as sample structures in a later chapter. 
However, the designer of an actual offshore structure, whose physical pro- 
perties and dimensions will eventually be known, must have a means of esti- 
mating the coefficients Cj and with reasonable accuracy in order to analyze 
the response of his structure. Further, investigators who are engaged in 
experimental research on the response of submerged or partially submerged 
bodies to earthquake excitation should have a method of estimating the coef- 
ficients in order that the experimental work may be conducted over a reason- 
able range of variables. Since Cj and are most commonly considered to be 
dependent upon the relative motion between the mass and the fluid, x $ and x $ 
in the case of the mass moving through the fluid at rest, a method is developed 
in this section for estimating the response quantities lx | , <|x |>, lx | , 

<|x s |> (hereinafter called hydrodynamic parameters) for strong motion earth- 
quakes of the type considered in this thesis. 

Consider the equations 

MZ + CZ + K n .Z + KZ = -Mx - K m x (2.6) 

S DL g DL g ' ' 



and 
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Z + 2$ o)Z + 2B k wZ + uj^Z = -x - 28 lO)X (2.7) 

s h g n g 

is which the nonlinear drag term Kg has been linearized in the manner dis- 
cussed previously. In the response computation, the quantities x g and x $ 
may be computed at each time point from the relations 



and 

x = Z + x 
s g 

Equation 2.7 is solved numerically for values of B g and 3^ which are reasonable 
to expect for the type of system considered, and again it is found that the 
hydrodynamic parameters l x s l max > < |x s |>, l x s l max » anc * < l x s l > are affected by 
the total amount of damping in the system and not by the distribution of damp- 
ing between structural and hydraulic. Comparisons were made of the hydro- 
dynamic parameters when the total amount of damping was equal to 2, 5, 7, 

10, 12 and 15 percent of critical with the damping being distributed between 
B s and 3^ as shown in Table 2.3. The difference in a hydrodynamic parameter 
at a given value of total damping ratio is within two to five percent for the 
different distributions of damping between B s and 3^, thus indicating that, 
if the amount of total damping can be estimated, the parameter can be esti- 
mated. The comparison holds true for all four hydrodynamic parameters. 

Spectra for the hydrodynamic parameters l x s l max > <|x s |>, l x s l max » 
and <|x s |> are shown in Figs. 2.27 through 2.30 for total damping ratios up 
to ten percent. The general shape of the spectra indicates that an idealized 
spectrum may be able to be constructed for a given set of conditions in much 
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the same manner as an idealized deformation response spectrum is constructed. 
Each hydrodynamic parameter spectrum was plotted, a smooth spectrum which 
approximated the actual spectrum was established by the author, and appro- 
priate amplification factors were computed for each of the four earthquakes 
considered. These amplification factors were then averaged in a mean square 
sense to yield the amplification factors listed in Table 2.5. The use of 
these amplification factors to construct idealized spectra for the hydro- 

dynamic parameters lx I „ , <|x |>, lx I, and <|x |> is explained in the 
J r 1 s 'max 1 s' 1 s 'max 's' 

next section of this chapter. 

2.7 Idealized Elastic Hydrodynamic Parameter Response Spectra 

The construction of idealized response spectra for the hydro- 
dynamic parameters discussed in Section 2.6 is very similar to the construc- 
tion of an idealized elastic deformation response spectrum. 

2.7.1 Displacement Parameter 

To construct an approximate response spectrum for I x s I max (° r < I x s I > ) 
for a given system: 

1. Determine the quantities K^, to, 3 g > and 3^ in the same manner 
as described in Section 2.5. 

2. Enter the appropriate section of Table 2.5 with the value of 
3 $ expressed as a percent and obtain the displacement ampli- 
fication factor. Multiply this amplification factor times x 

r J r g max 

(or <|Xg|>) and plot the result as a displacement line on tri- 
partite logarithmic paper. This line is shown as line 1 on 
Fig. 2.31. 
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3. Enter the appropriate section of Table 2.5 with the amount of 

total damping (8 S + 8^) and obtain the velocity amplification 

factor. Multiply this amplification factor times x„ „„ (or 
r r g max 

<|Xg|>) and plot the result as a velocity line as shown by 
line 2 on Fig. 2.31 . 

4. Plot the value x „„ (or <|x |>) as a displacement line. This 

g max 1 g 1 

is shown as line 3 on Fig. 2.31. 

5. The envelope formed by the intersection of these three lines 
is an idealized spectrum for I x s l max (° r < l x s l > )- 

2.7.2 Velocity Parameter 

To construct an idealized response spectrum for l x s l max (or <|x s |>) 
for a given system: 

1. Obtain and use displacement and velocity amplification factors 
from the appropriate section of Table 2.5 exactly as described 
above for the construction of an idealized spectrum for a dis- 
placement hydrodynamic parameter. The resulting lines are shown 
as lines 1 and 2 on Fig. 2.32. 

2. Enter the appropriate section of Table 2.5 with the amount of 
total damping (8 S + 8^) and obtain the acceleration amplifica- 
tion factor. Multiply this amplification factor times x 

g max 

(or < | x g | > ) and plot the result as an acceleration line. This 
is shown as line 3 on Fig. 2.32. 

Plot the value of x^ max (or < | Xg | > ) as a velocity line as shown 
by line 4 on Fig. 2.32. 



3 . 
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4. The envelope formed by the intersection of the four lines is 

an idealized response spectrum for |x s | mgx (or <|x s |>). 

Linear interpolation can be used to obtain amplification factors 

for values of damping not shown exactly in the table. 

Unfortunately, no reasonably consistent ratios were found to exist 

between < | x^ | > , <|x g |>, < |x g |> and other ground motion quantities such as 

x , x „ , x _ for the four earthquakes considered. Thus, the de- 
g max g max g max ^ 

signer of an offshore structure who wishes to estimate the hydrodynamic param- 
eters <|x |> and <|x I > will have to obtain the values of<|x |>, <|x |>, and 
's' 's' 1 g 1 1 g 1 

< |Xg| > from processing an actual earthquake record or estimating them using 

the values shown in Table 1.1 as a guide. The range of the ratios between 

<|xl>, <|xl>, <|xl> and x„ m , x„ m . x M mav are repeated here for the 
1 g 1 1 g 1 1 g' g max g max g max 

convenience of the reader. 



Ratio 



x g max 
x g max 
x g max 



Lowest Value 
0.194 

0.124 

0.084 



Highest Value 
0.385 

0.167 

0.185 



It should be noted that the idealized spectra for hydrodynamic param- 
eters become less descriptive of the oscillator response at very low frequencies, 
about 0.08 cps and less for the four earthquakes used in the computations. This 
is to be expected, since the quantities lx | , <|x |>, |x | , and <|x |> 

become quite small for very soft (low frequency) systems. 
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2.8 Inelastic Response Spectra 

In this section, the effects of the nonlinear hydrodynamic drag 
force on the response of the oscillator in the inelastic range are investi- 
gated. Figure 2.33 shows the resistance-deformation relationship of the 
spring in a system with elasto-plastic spring characteristics. The character- 
istics of the elasto-plastic system are an initial elastic region, a constant 
plastic resistance limit, and an elastic unloading region. The unloading 
portion is considered to have the same slope as the initial elastic curve, and 
the unloading is elastic until yielding is reached in the opposite direction. 12 
The yield point deformation is denoted by Z^, and the absolute maximum dis- 
placement of the mass relative to the base is denoted by | Z l max * The yield 
resistance is denoted by R . The ductility factor, y, is defined as 




and expresses the maximum deformation of the inelastic system in terms of its 
yield point deformation. With y defined in this manner, the maximum deforma- 
tion of the system is yZ^, the maximum elastic deformation is Z y , and the 
maximum inelastic deformation is (y - 1 )Z y . A variety of calculations of 
the response of inelastic systems has been made and reported in the litera- 
ture. 11-14 The results of the calculations are conveniently displayed in 
the form of response spectra in much the same manner as response spectra 
for elastic systems are displayed. For elasto-plastic systems, response 
spectra can be plotted for constant values of the ductility factor, y, and 
when the plot is made in terms of the elastic component of the displacement, 
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the accelerations are properly presented. 12 In this type of spectrum represen- 
tation, the pseudo-velocity quantity plotted is u)Z y . Elasto-plastic response 

spectra for systems with 6 = 0.02, Lx n m / K = 0.00 (no drag), subjected 

s u g niaX 

to the El Centro earthquake are shown in Fig. 2.34 for values of p = 1 (elastic 
case) to y = 4. In these spectra, the total deformation of the system can be 
determined by reading the displacement quantity, which is Z y , and then using 
the relationship 

I ^ 1 max ^ y 

In the low frequency range of the spectra shown in Fig. 2.34, the 
total deformation is approximately the same for the elasto-plastic systems 
as for the elastic systems. In the high frequency region, the pseudo-acceler- 
ation is approximately the same for the elasto-plastic systems as for the 
elastic system. In the mid-frequency region, some generalizations have been 
made regarding the amount of reduction in yield level for different values 
of y, but more detailed generalizations are still under study. 12 

Elasto-plastic deformation response spectra for systems with non- 
linear hydrodynamic drag forces subjected to the El Centro earthquake are 
shown in Figs. 2.35 through 2.37. Again, in the low frequency region, the 
total deformation is approximately the same for the elasto-plastic systems 
as for the elastic system. In the mid-frequency region, the spectra have 
the same general shape as do the elasto-plastic spectra for systems without 
the drag forces, but no quantification of the results has been made. Further 
study in the direction of establishing rules for the construction of idealized 
elasto-plastic response spectra in this frequency region is clearly indicated. 
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In the high frequency region, the elasto-plastic spectra for the systems with 
the lowest value of the drag parameter (0.006) approach approximately the same 
values in the high frequency tail as the elastic spectrum does. However, for 
the systems with the greatest value of the drag parameter (0.50), the high 
frequency tails of the elasto-plastic spectra approach values which are re- 
duced by as much as 50 percent from the elastic case. In other words, there 
is a reduction in the magnitude of the elasto-plastic high frequency tail from 
that obtained in the elastic case, and the amount of the reduction is a func- 
tion of the amount of nonlinear hydrodynamic drag in the system. In the sys- 
tems with a small amount of drag, the reduction is very slight, and may not 
merit much consideration, but the reduction for high drag systems can be sig- 
nificant. A study of the numerical results obtained in this investigation 
revealed that, for the spectra with 8 = 0.02, K n x ,/K = 0.50, values of 

p up to 4 and the El Centro earthquake, the factor y applied to the 
elastic case high frequency tail yields a satisfactory estimate of the elasto- 
plastic high frequency tail (see Table 2.6). The reader is cautioned that 
this expression for the reduction factor is applicable only to the specific 
systems and earthquake used, and is not to be extended to other conditions. 

In order to compare the effects of the nonlinear hydrodynamic drag 
force on the elasto-plastic response spectra, spectra are plotted for given 
values of structural damping ratio, 6 S » and ductility factor, p, for systems 
with varying amounts of drag subjected to the El Centro earthquake. These 
comparisons are shown in Figs. 2.38 through 2.40. These figures show that, 
for a small ductility factor (p = 1.5), increasing amounts of nonlinear hydro- 
dynamic drag reduce the value of the yield deformation, Z^, for a range of 
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frequencies between low frequencies and the beginning of the high frequency 
tail. This means that the maximum response of the system with high drag forces 
is less than the maximum response of the system with lower drag forces, or, 
the drag forces tend to produce a damping effect on the maximum deformation 
of the spring in the mid-frequency region of the spectrum. For example, in 
Fig. 2.38, at a frequency of 1 cps, and the highest value of the dimensionless 
drag parameter used (0.50), is only about 50 percent of the value of 
for the system with no drag. However, as the ductility factor y increases, 
the effects of the nonlinear hydrodynamic drag force become less and less 
noticeable in the mid-frequency range. In Fig. 2.40, with y = 4, the value 
of Zy for the highest drag system is about 90 percent of the value of Z^ for 
the system with no drag at a frequency of 1 cps. This is reasonable to expect, 
since the velocity of the oscillator, x s , can be expected to be smaller during 
plastic deformation of the spring than during elastic deformation of the 
spring, and the magnitude of the drag force is a function of the oscillator 
velocity. It is also to be noted from Fig. 2.40 that it is possible for Z^ 
in the case of the system with drag to exceed slightly the value of Z y for 

the system with no drag at the same frequency (note the response values at 

frequencies of 0.8 cps and 3.0 cps) in the mid-frequency range. This can re- 

sult since the maximum response of the elasto-plastic system depends on the 
entire history of the deformation and the energy absorbed in the spring during 
plastic deformation as well as the energy dissipated by the damping effects 
of the fluid surrounding the oscillator. 

The significant features of the elasto-plastic response spectra 
for systems with nonlinear hydrodynamic drag are: 
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1. The drag does not have much effect on the spectrum in the low 
frequency region; 

2. A high frequency tail does exist, and large amounts of drag 
can cause significant reductions in Z for increasing values 

y 

of the ductility factor (less significant reductions in with 
increasing values of y for low drag systems); 

3. The damping effects of the drag in the mid-frequency region 

of the spectrum become less pronounced as the ductility factor 
increases, even for the systems with high drag. 
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CHAPTER 3 

MULTI-DEGREE OF FREEDOM SYSTEMS 



3.1 General 

In this chapter, the principles developed and discussed in Chapter 
2 for the single degree of freedom oscillator are extended to the multi- 
degree of freedom system. The nonlinear drag forces are linearized in such 
a manner as to insure the existence of classical normal modes in the equa- 
tions of motion in order that normal mode superposition may be used as an 
approximate solution to the problem. The method of linearization is shown 
to be accurate for a series of five degree of freedom systems which are 
representative of the type of structure being investigated. A means of es- 
timating the amount of equivalent linear hydraulic damping, using ground 
motion velocity characteristics, is developed and is shown to yield good re- 
sults. The construction of an idealized response spectrum for each mode of 
the multi-degree of freedom system is discussed and described, and sample 
problems are solved by the superposition of normal mode responses. 

3.2 Approximate Solution of the Nonlinear Equations of Motion 

In this section, an approximate method of solving the nonlinear 
equations of motion for the multi-degree of freedom system is developed and 
shown to yield satisfactory results. Various features of the solution, such 
as the formulation of the structural damping matrix, the linearization of 
the drag forces, the superposition of modal responses, and the response of 
very high frequency modes, are discussed. The transformation of the nonlinear 
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equations of motion into uncoupled linear modal equations, which was first 
shown in Chapter 1, is repeated in this section for the convenience of the 
reader. 

3.2.1 Equations of Motion 

The equations of motion for the lumped parameter model shown in 
Fig. 1.1(b) are written in matrix form as Eq. 1.3 and are repeated here to 
form the basis of the analytical discussion in this section. 

C-M-J{Z> + [C s ]{Z> + [K]{Z> = -FMj{x g ) - rK D J (x s |x s | } (3.1) 

No structural dampers are shown on the model, but the damping matrix [C s ] 
can take a variety of forms as will be discussed in later paragraphs. For 
a given structural system and known earthquake record, response quantities 
of interest may be obtained by the numerical integration of Eq. 3.1, but 
this procedure is costly and requires the use of a digital computer. If 
Eq. 3.1 can be satisfactorily approximated by a system of simultaneous 
linear differential equations of motion which possess classical normal 
modes, then a solution can be obtained by the superposition of normal mode 
solutions. A further approximation can be made by using response spectra 
to obtain the maximum response in each mode and then combining the modal 
maxima in some fashion, very much like the approximate modal techniques, 
using response spectra, which are reported in the literature 11-13 for land- 
based structures. 

In the normal mode method, the vector of a response quantity ex- 
pressed in one coordinate system, such as the relative displacement between 
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each node and the base, {Z}, is transformed into the "normal" coordinate sys- 
tem by the relation 

{Z> = [4>]{Y> (3.2) 

If the nonlinear drag force were not present, Eq. 3.2 and the ortho 
gonality relationship between [<f>] and ["M.J, [K], and [C s ] could be used to 
transform Eq. 3.1 into nm (number of modes) uncoupled single degree of free- 
dom equations in the normal coordinate system. If the nonlinear drag force 
can be converted into an equivalent linear drag force, then the resulting 
linear differential equations of motion can be transformed into uncoupled 
single degree of freedom equations of motion in the normal coordinate system 
as follows: Consider that C'Kp-] {x s |x 1} can be transformed into an equiva 

lent linear force, {x $ } i.e., 

[■'K D J {x sl. x s |} - n< DL J {x $ } (3.3) 

Substitution of Eq. 3.3 into Eq. 3.1 yields, 

rMj{Z} + [C $ ]{Z}+ [K]{Z> = -C'Mj{x g } - IT~K DL -J{x s } (3.4) 

Since {Z} = {x s > - {Xg> , Eq. 3.4 can be rewritten as 

rM J{Z} + [C S ]{Z> + rK DL J{Z} + [ K] { Z} 

= -FMjfXg} - rK DL j{x g > 

Substitute Eq. 3.2 into Eq. 3.5 and premultiply the result by [4>]^ 



(3.5) 
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M T t-MJM{Y} + W T [c s ][*](n * mYk dl JM{y} + W T [K]W(Y} 

■ - M T rK DL J<x g > (3.6) 



Recall the following definitions from Chapter 1: 



* 

D-HJ 



* 

CK J 






o 2 j 



M T MM 

[*] t [c s ]M * 

M t [k]M 
m t c-k dl jm = 



a diagonal modal mass matrix 

a diagonal modal structural damping 
matrix. The form of [C s ] which permits 
this is discussed later, 
a diagonal modal stiffness matrix 

a modal linear hydraulic damping 
matrix which is not necessarily 
diagonal 

a matrix of circular natural fre- 
quencies of the system 



Substitution of these definitions into Eq. 3.6 yields 

['M*J{Y> + + [Kp L ]{Y} + [" K*~J { Y } 

= -fo] T FMj{X g } - [4»] T rK DL J{x g } (3.7) 

At this point, the formulation of the structural damping matrix, [C ], and 
the linearization of the drag factor, will be discussed. 

3.2.2 Structural Damping 

If the structural damping matrix [C ] is proportional to the mass 
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matrix, then [“C^] will be diagonal since the natural mode shapes are ortho- 
gonal with the mass matrix. This formulation of [C s ] results in modal damping 
ratios which decrease with increasing natural frequency. This is shown as 
follows: 

Let, 

[C s ] = SFMJ 



then, 

m t [c s ]m = sofr-M-jw = = c-c*j 

Define 

= r2e s wj = £[I] 

then, 

8 = _L_ 

s i 2 “i1 

where 

3 = structural damping ratio in the ith mode 

s i 

= circular natural frequency of the ith mode 



( 3 . 8 ) 



( 3 . 9 ) 



( 3 . 10 ) 



In a similar manner, it can be shown that, if the structural damping 

‘A’ 

matrix [C ] is proportional to the stiffness matrix [K], then [~ C is a 

diagonal matrix in which the modal damping ratios 3c increase with increasing 

* i 

frequency. 

The elements of the structural damping matrix [C s ], can be selected 
to yield the same structural damping ratio in each mode. It is common prac- 
tice to assume the same structural damping ratio in all modes in the dynamic 
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analysis of damped multi-degree of freedom systems. 12 If the modal matrix 
[4>] is normalized such that 



[<f>]^[' M ~J[<f>] = F I J (a unit matrix) 



(3.11) 



Then 



[C s ] = rMJ[4,]rC*J[4,] T C-MJ 



(3.12) 



And, since 



= r'M^IT 2B s 0i-J 



(3.13) 



Equation 3.12 can be written as 



[c s ] = rM j[<f>]rM*^ir26 s mj[<f)] T rMj 



(3.14) 



Thus, it is seen that [C ] can be formulated in such a manner as 
to yield modal structural damping ratios which are constant, which increase 
with increasing mode number, or which decrease with increasing mode number. 

In the computations reported in this thesis, the structural damping matrix 
[C ] is computed in such a way that is the same in all modes. 

3.2.3 Linearization of the Drag Force 

Foster, 9 and Malhotra and Penzien 9 linearize the nonlinear drag 
forces, [■'Kp-Jtx |x |} by replacing them with a velocity dependent linear 
force, ["K|^~~]{x s } at each node which performs the same amount of work during 
a period of time as does the nonlinear force at that node. Foster then 
recognizes that the resulting matrix is not orthogonal with the mode 

shapes and proceeds to solve the equations of motion by a nonclassical mode 



57 



approach using complex eigenvalues and eigenvectors. Malhotra and Penzien 

^ r\j 

transform into [K^] and then force [K^] to be diagonal by means 

of an iterative process that replaces the nondiagonal matrix with a diagonal 
matrix which dissipates the same amount of energy over a period of time. 

The hydraulic drag forces, whether linear or nonlinear, may be con- 
sidered to act as dampers, that is, they cause energy to be dissipated as 
the system oscillates through the fluid. In the approximate solution being 
developed in this section, the nonlinear dampers in the multi -degree of 
freedom system are replaced with linear dampers (except for very high fre- 
quency modes which are discussed in a later section) which dissipate the 
same amount of energy over the total system, but do not necessarily preserve 
the dissipation of energy at each node. Further, the resulting "equivalent" 
linear damper at each node is made to be proportional to the lumped mass of 
that node in order that the modal equations of motion will be uncoupled. For 
an offshore structure with a portion of the structure above the water line 
as shown in Fig. 1.1(a), this procedure results in a linear hydraulic damper 
being provided at the level of the top mass where, in fact, no nonlinear 
hydraulic damper really exists. With this exception, the proportioning of 
the linear hydraulic dampers in the same ratio as the mass is distributed is 
a reasonable approach for the type of structure being discussed, since the 
nonlinear damper is proportional to the projected area of the structure lumped 
at the node, and the magnitude of the mass at a given node will also be gen- 
erally proportional to the projected area. The method of linearizing is as 
follows: 
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Energy dissipated by 
nonlinear dampers 



l 

j=l 



(k d.. vivi ){ v } dt 
jj j j j 



= t i 
j=i 





(3.15) 



Energy dissipated by nn /-t . 

equivalent linear = £ (1C. x )(x ) dt 

dampers j=l £ UL jj s j s j 



nn ? 

- 1 l K DL ..<> < 3J6 > 

J*1 JJ J 



where 



nn = number of nodes in the system 
j = jth node in the system 



Equating the energy dissipated. 



nn 



l 

j=l 








(3.17) 



It is noted that only the total energy dissipated is preserved and the dis- 
tribution of the energy dissipation among nodes is not preserved. Now, by 
forcing the linear dampers to be proportional to the mass matrix. 






(3.18) 



and substituting Eq. 3.18 into Eq. 3.17 
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a 



nn .? . 

.1 K D.. <x sJ x s.I > 
J=1 JJ J J 



nn 

I M, 



<x 2 > 
j=l JJ s j 



(3.19) 



% T 

Recall that [K^] was defined as [<}>] C" K DL —J [4>] . Substituting Eq. 3.18 into 
this definition. 



[K dl ] = otdd)] 1 ^ M = a^M 



(3.20) 



Substituting Eq. 3.20 into Eq. 3.7 yields 

rM*^]{Y} + d'C*J{Y} + aC'M*J(Y} + |M<\]{Y} 

= C4>] T r M J{-x g - ax g } (3.21) 

Equation 3.21 represents nm single degree of freedom uncoupled equations of 
motion in the normal coordinate system. 

Before proceeding, the method of linearizing described by Eqs. 3.18 
and 3.19 will be tested for accuracy in computing the response of multi -degree 
of freedom systems to earthquake excitation. A series of five degree of free- 
dom systems, representative of offshore towers in mass, stiffness, and drag 
distribution, is developed in order to test the linearization technique. 

Seven structural systems, consisting of different combinations of 
mass and stiffness distribution, are considered. Each system is further 
subdivided into one or more distributions of the drag factor, K^, at each node. 
The characteristics of these systems are shown in Table 3.1. In the table, 
the values shown for the mass at each node are considered to be the total 



mass at that node, and include the added mass term. The different 
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distributions of drag factor are labeled L, M, and H in the table, and the 
column under each of these headings represents the relative magnitude of 
at each node for the particular distribution shown. For example, drag factor 
distribution M represents the distribution of at each node as 0, 1, 1, 2, 

3 at nodes 1 through 5 respectively. The values of ZKp/Zm shown in the 
table for the distributions of drag factor are of sufficient range to span 
reasonable values of this quantity for actual offshore towers of the type 
being investigated. The quantity EK^/Em for the four towers reported by 
Burke and Tighe 4 range from a low of 0.0034 inches to a high of 0.0054 
inches - ^. The relative magnitudes of the undamped natural frequencies of 
the system, normalized to the fundamental frequency, are as shown. In the 
computations which are discussed in the following paragraphs, one of the 
natural frequencies is set to a desired value and the fundamental frequency 
is computed from the ratio of natural frequencies. The value of is com- 
puted from f.j and the quantity f-|/A-|/nii , then the remaining spring stiffnesses 
are computed from k-j and the distribution of k ^ through k^. The stiffness 
matrix [K] for the entire system is then assembled and is a tri-diagonal 
matrix for the systems shown. 

Each of the systems shown in Table 3.1 is subjected to the El Centro 
earthquake and response quantities are computed by the numerical integration 
of the equations of motion using the criteria for convergence and time inter- 
val described in Chapter 1. For each system, computations of maximum re- 
sponse are made for two values of in each mode, 0.01 and 0.05, and three 
frequency ranges. The low frequency range, range number 1, is established 
by setting the fifth natural frequency equal to the frequency at which x^ mQx 
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a nd Xg max intersect when plotted on tripartite logarithmic paper. The 
mid-frequency range, range number 2, is established by setting the funda- 
mental frequency equal to the frequency at which x and x„ mziv inter- 

9 iTlaX 9 

sect. The high frequency range, range number 3, is established by setting 
the fundamental frequency equal to the frequency at which x^ max and x^ max 
intersect when plotted on tripartite logarithmic paper. The maximum dis- 
placement of each node relative to the base, {|Z| }, is computed. The 

results of the computations are shown in Tables 3.2 through 3.27. The 
column in the tables titled "No Drag" shows the results of the computations 
when ["KpJ is set to zero. The column in the tables titled "Actual Drag" 
shows the results of the computations when ["Kp~J is computed from the values 
shown in Table 3.1, thus nonlinear hydrodynamic drag forces, [~ K^-JtxJxJ }, 
are acting on the structure. In the computations for the "Actual Drag" 
systems, the time averages {<x^|x s |>} and {<x^>} are collected. These time 
averages are then used to linearize the drag factor and set it to be pro- 
portional to the mass matrix as prescribed by Eqs. 3.18 and 3.19. Equation 
3.4, using the linear drag factor matrix just computed, is then integrated 
numerically to obtain the response quantity {|Z| max K The results of this 
integration are labeled "Case 1" in the tables of results. The column 
labeled "Case 1/Actual Drag" in Tables 3.2 through 3.27 represents the com- 
parison of the solution of Eq. 3.1 with the solution of Eq. 3.4. A study 
of these results reveals that the linearization technique described by Eqs. 
3.18 and 3.19 is extremely accurate. The difference between the linear and 
the nonlinear results is less than 6 percent (with one exception — system 1M, 
B s = 0.01, node 3, mid-frequency range) and is less than 2 percent in the 
majority of the cases. 
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The maximum distortion between two adjacent nodes is also computed 
for both the linear and the nonlinear equations of motion. The results of 
the computations for the linear case in which the linearization is carried 
out according to Eqs 3.18 and 3.19 are not shown in this thesis, but the 
linear solution is as close to the nonlinear solution as for the response 
quantity {|Z| max } shown in the tables. Thus, it is concluded that the meth- 
od of linearizing the drag force described by Eqs. 3.18 and 3.19 yields 
satisfactory results. 

The primary difficulty in using Eqs. 3.18 and 3.19 to linearize the 

drag force lies in the ability to determine the quantities {<x^|x s |>} and 

• 2 

{<x s >} a priori . If these quantities can be estimated in some fashion, and 
if the estimate can be shown to be accurate in certain frequency ranges, 
then Eq. 3.1 can be linearized into Eq. 3.4, and Eq. 3.4 can be transformed 
into uncoupled, modal, single-degree of freedom equations in the normal co- 
ordinate system which can be solved by spectrum techniques. It was demon- 
strated in Chapter 2 that, for the single degree of freedom oscillator, the 
quantities <xf|x l> and <x > are conservative estimates of the oscillator 

g 1 g 1 g 

velocity in all but the low frequency region of the spectrum. Since the 
major effects of hydraulic damping occur in the velocity amplification re- 
gion of the spectrum, it seems reasonable to estimate the quantities 

* 2 * *2 * 2 * *2 

(<x | x | >} and {<x^>} as {<xnx |>} and {<x >} respectively. Using these 
s s s 9 9 9 

quantities, Eq. 3.19 becomes 



/ nn 

- 1 , K| 

J=1 



D. . 
JJ 



^ / • 2 1 • | \ 

<x g |x gl > 



nn 

l M. . , 

' j=l ^ / 



<x g 



a 



(3.22) 
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A series of computations, labeled "Case 2" in Tables 3.2 through 3.27, is 
made in which the nonlinear drag factor is linearized in the manner pre- 
scribed by Eqs. 3.18 and 3.22. A study of the column labeled "Case 2/Actual 
Drag" shows that the linearization described by Eqs. 3.18 and 3.22 is, in 
general, conservative in the mid and high frequency ranges and is generally 
unconservative in the low frequency range. The degree of conservatism in 
the mid-frequency range is seen to be less than 5 percent in most cases, 
although one case is conservative by 12 percent. In the high frequency 
range, the results are seen to be conservative by as much as 20 percent in 
some cases, especially in the high drag systems, but many of the results are 
conservative by less than 10 percent in the high frequency range. The degree 
of unconservatism in the low frequency range is seen to be less than 10 per- 
cent in most cases, but can range as high as 23 percent. This unconservatism 
in the low frequency range is to be expected, since very soft (low frequency) 
systems can be expected to exhibit little velocity response at each node, 
and so the use of the ground velocity characteristics overestimates the 
velocity of each node. One further observation to be made from a study of 
the numerical results is that the use of the ground velocity to linearize 
the drag force is generally more conservative as the response is reduced 
to a greater degree by the actual nonlinear drag force. In those cases in 
which the Case 2 linearization is unconservative in frequency ranges 2 and 
3, the degree of unconservatism is generally very slight (less than 1 percent). 

The maximum deformation between adjacent nodes is also computed 
and a comparison of the linear drag response (linearized by the ground velocity 
characteristics) with the nonlinear drag response is shown in Figs. 3.1 
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through 3.13. These figures show that the maximum spring deformation com- 
puted when the drag force is linearized by the ground velocity characteristics 
is also conservative in the mid and high frequency ranges and generally un- 
conservative in the low frequency range. It is concluded from these compu- 
tations that the use of Eqs 3.18 and 3.22 is a reasonable way to linearize 
the drag force and yields slightly conservative results except in the low 
frequency range. The conservatism introduced by the linearization is seen 
from the figures to be generally less than 10 percent in the mid and high 

frequency ranges, with some of the results in the high frequency range being 

conservative by as much as 20 percent. The results of the low frequency range 
calculations are generally unconservative by less than 20 percent, but can be 
unconservative by as much as 25 percent. 

The effect of structural damping on the response is quite apparent 
from the tables and figures. For example, the maximum relative displacement 
between the bottom node (mass number 5) and the ground for system 4M in fre- 
quency range 2 is reduced approximately 17 percent by the hydraulic drag when 
the structural damping ratio is 0.01 in each mode (see Table 3.14). Yet the 
response is reduced by less than 2 percent by the hydraulic drag when the 

structural damping ratio is 0.05 in each mode (see Table 3.15). Thus, the 

selection of structural damping ratios to use in the response computations 
is seen to be very important. 

3.2.4 Modal Equations in Normal Coordinates 

The use of Eqs. 3.18 and 3.22 to linearize the drag force has been 
shown to yield generally conservative results in the mid and high frequency 
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ranges, and unconservative results in the low frequency range. Since the 
linearization prescribed by Eqs. 3.18 and 3.22 insures that the equations of 
motion for the linearized system (Eq. 3.5) will possess classical normal modes. 



the discussion of the superposition of normal modes is resumed, 
uncoupled equations of motion in the normal coordinate system. 

C'M*^KY} + FC*J{Y} + a["M*J{Y} + ^K*J{Y} 


Recall the 


= [*] T FMj{-x g - ax g } 

ic _ I 

Premultiplying by ['M yields 


(3.21) 


m + [-26 s wjm + [~2e h u.-]m + o 2 J{y> 

= rH*-J' 1 W T CMj{-x g - ax g } 

where 


(3.23) 


ari-j = r 23 h a>J 

The terms on the right side of Eq. 3.23 can be written as 


(3.24) 


r'M*^]" 1 [4.] T ['M J(-x g - otXg} = {PF} (-x g - aXg) 


(3.25) 


where 

{PF} = ['M*^ _1 [4»] T r'Mj{I} 


(3.26) 



and PF.. is the participation factor for the ith mode. Substituting Eq. 3.26 



into Eq. 3.23 

{Y} + r26 s wJ{Y} + T26 h a)J{Y} + |> 2 J{Y} 
= {PF}(-x g - ax g ) 


(3.27) 



(3.27) 
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For the ith mode, then, 

Y, ♦ ZB^Y, + + ^Y, ■ "V'*, - «,) < 3 ' 28 » 

Since the left side of Eq. 3.28 is linear with respect to the unknown Y, the 
maximum modal response in the ith mode can be obtained by determining the 
maximum response of the single degree of freedom equation. 



y. + 28 a), .y . + 28 l. w-*y. + to..y. 
p s i lri p h. lri iri 



-x - ax 

g g 



(3.29) 



and multiplying this response times the participation factor for the ith mode, 
PF. . The nonlinear equations of motion for the multi-degree of freedom system, 
Eq. 3.1, have now been transformed into a set of linear simultaneous differ- 
ential equations (Eq. 3.5) by the use of Eqs. 3.18 and 3.22. The linear 
simultaneous differential equations are transformed into nm uncoupled, single 
degree of freedom modal equations of the form represented by Eq. 3.28 in the 
normal coordinate system Y by the use of Eq. 3.2. The total response of the 
system is the result of a combination of the various modal responses represent- 
ed by Eq. 3.28. 



3.2.5 High Frequency Modes 

The existence of a high frequency tail in the response spectrum 
for a single degree of freedom oscillator was demonstrated and discussed in 
Chapter 2. Further, the linearization of the nonlinear hydrodynamic drag 
force was shown to break down when this high frequency tail began to appear. 
A method for determining the parameters necessary to estimate the magnitude 
of this tail for each mode in the multi -degree of freedom system is now 
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developed. Consider the equations of motion for an undamped multi-degree of 
freedom system 

rMjm + [k]{z> = -['Mj{x g > - r-K D j{x s |x s |} 0.30) 

If the system is so stiff that {x } ~ {x }, Eq. 3.30 can be written as 

s g 

C'MjtZ) + [K]{Z} = -C" M -J{x g ) - rK D J {x g |x g | } (3.31) 

Transform Eq. 3.31 into the normal coordinate system by the use of Eq. 3.2 
and premultiply by [<j>] T 

W T fMO[*]{V) + [*] T [K]M{Y) 

* -[*] T ['MjtI}x g - [*] T C'K D J{I}(x g |J g |) (3.32) 

or, 

+ C'K*J{Y} 

■ -[♦] T ['Mj{I}5i g - W T CK D J{I)(J g |J g |) (3.33) 

Equation 3.33 represents nm single degree of freedom uncoupled modal equations 
in normal coordinates for the high frequency multi-degree of freedom system. 
Premultiply Eq. 3.33 by then 

m + r<Ajm = -{PF>x g - {y(£ g |x g p (3.34) 

where 

{K*} = rM*J' ] [())] T rK D J{I} (3.35) 



For the ith mode then. 
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K 



D. 



Y. + o)..Y. = PF. (-x - DC /\ i^i 

i n l l g PF.. g 1 g 1 



^xjxl) 



(3.36) 



Again, since the left side of Eq. 3.36 is linear with respect to the unknown 
Y, the maximum modal response in the ith mode can be obtained by determining 
the maximum response of 



.. . d 

y . + a), .y . = -x 

J -\ lri c 



PF. 



V V 



(3.37) 



and multiplying this response times the participation factor for the ith mode. 
It was shown in Chapter 2 that the quantity |y . | can be estimated for 

I IlldX 

Eq. 3.37 as 



l y i 



max 



K, 



°i 



\ 



PF. 



9 * 

o)T .M. . 
\ n n 



( *g max*' 



(3.38) 



and 



i 'max 



’ < PF i><l*il m ax> 



(3.39) 



3.2.6 Superposition of Modal Responses 

The maximum response in any given normal mode of the multi -degree 
of freedom system can be computed from Eq. 3.28. However, the maximum re- 
sponse in one mode does not necessarily occur at the same time as the maximum 
response in another mode, and the time of occurrence of the maximum response 
in the time history of the modal response is not normally known. The actual 
combination of modal responses in a multi -degree of freedom system is very 
complex and few generalities can be made about the combination. 
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An approximation to the actual maximum response of a multi-degree 
of freedom structure can be made by determining the maximum response in each 
mode, and then combining the maximum modal responses in some fashion. An 
upper bound on the maximum response of the structure can be established by 
assuming that the maximum modal responses occur at the same time in each 
mode and are additive in nature. In other words, the summation of the ab- 
solute values of the maximum response in each mode yields an upper bound to 
the response of the structure. Another commonly used method of combining the 
maximum modal responses is to take the square root of the sum of the squares 
of the modal responses. Other methods of combining modal responses are some- 
times used. 22 The two most commonly used methods are the summation of the 
absolute values of the modal maxima (SABS) and the square root of the sum 
of the squares of the modal maxima (SRSS). These two methods are the only 
ones used in this thesis for illustrative purposes. Tables 3.28 through 
3.33 show the results of the superposition of maximum modal responses applied 
to three of the five degree of freedom systems. The drag force is linearized 
in accordance with Eqs. 3.18 and 3.22 (ground velocity linearization), and 
the amount of linear hydraulic drag B h in each mode is computed from Eq. 3.24. 
The maximum modal response in each mode is computed from the response of the 
single degree of freedom modal equation (Eq. 3.29), and multiplied times the 
participation factor for the mode. The EABS combination is seen to be higher 
than the true response (as should be the case) and the SRSS combination is 
generally less than the true response. No generalizations as to the relative 
accuracy for different distributions of mass, stiffness, and drag factor can 
be made, however. 13 Thus, it is seen that, if the maximum responses of the 
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single degree of freedom modal equations, Eq. 3.28, can be computed, the maxi- 
mum response of the structure can be approximated by a combination of the 
modal maxima. 



3.3 Idealized Elastic Deformation Response Spectra 



It was seen in the preceding section that the maximum response of 
a multi-degree of freedom system can be approximated by combining maximum 
modal responses in some fashion. If the maximum modal responses can be esti- 
mated in a way that accounts for the effects of the nonlinear hydrodynamic 
drag forces in the system, a rapid approximation of the response of the sys- 
tem can be obtained. The principles developed in Chapter 2 for the single 
degree of freedom oscillator and the method of linearizing the drag forces 
developed in this chapter can be combined to construct an idealized elastic 
deformation response spectrum for each mode of the multi -degree of freedom 
system subjected to strong motion earthquake excitation. The procedure for 
the construction of the modal elastic deformation response spectra is now set 
forth. 



1. Compute the natural frequencies, ["wj, and mode shapes, [<})], 
for the system. 

2. Establish the structural damping ratio, B s ., for each mode. 
Values of Be. between 0.005 and 0.05 are reasonable to assume 
for the type of structure being investigated in this thesis. 

3. Determine x , x , x , and 

g max g max g max 



Slid: 

4 
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for the earthquake to be used. If the spectra are to be con- 
structed for a hypothetical or "design" earthquake, one-half 

of the value of x is a reasonable value to assume for the 
g max 

quantity 



4. 



i\,i> 

■k k 

Compute M^. , Kp , PF.. , a, 3^., and the high frequency tail 
estimate from the following expressions: 



t~M*j = mVmjm 



(U - C'M*J‘ , W T C'K n JU} 



{pf} = jm 



(3.35) 

(3.26) 



a = 



/ nn 1 
.1 s.. 

J=1 JJ 




•2,. | \ 
<X 9 |X 9 I> 
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I M.. i 




X 9 / 


\ j=l JJ / 




/ 



(3.22) 



r26 h co-J = a [Mj 



(3.24) 



High Frequency Tail 
Estimate for the ith = 
Mode 



,n D. 

i 

PF. 



9 k 

“ii M n 



<*g 



(3.38) 



With the quantities to. . , 3 S .» 3h , and the high frequency tail 

11 i i 

estimate, construct the idealized elastic deformation response 
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spectrum for the ith mode exactly as described in Section 2.5 
for the single degree of freedom system. 

It is noted that the idealized response spectrum constructed as 
outlined above accounts for the fluid-structure interaction in a slightly 
conservative manner. When an idealized spectrum has been drawn for each mode 
of the system, the maximum modal responses are obtained from the spectra, 
multiplied times the appropriate modal participation factor, and combined in 
some manner to yield an approximate maximum response of the system. 

3.4 Hydrodynamic Parameters 

An approximate method for computing hydrodynamic parameters for a 
multi -degree of freedom system by the superposition of normal mode responses 
is developed in this section. As explained in Chapter 2, these hydrodynamic 
parameters are absolute displacement and velocity responses of the mass (each 
node in the multi-degree of freedom system), whose magnitudes are desirable 
to know in order that the empirical drag and inertia coefficients, Cq and Cj, 
may be estimated in a reasonable way. The hydrodynamic parameter spectrum 
amplification factors developed in Chapter 2 may be used to construct idealized 
spectra from which the response quantity of interest in each mode can be esti- 
mated. 

Consider again the linear equations of motion, 

C'MJ{Z} + [C s ]{Z> + [ K] { Z} = -["M-v]{x g } - rK DL J{x s ) (3.4) 

Making the substitution {Z} = {x s > - {x g }, Eq- 3.4 can be written as 
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+ [C $ ](x s } + rK DL J{x s } + [K]{x s > 

= [K]{x g } + [C s ](x g } (3.40) 

If the velocity force, [C s ]{x g ], can be dropped from Eq. 3.40 for values 
of [C ] up to some magnitude without unreasonable loss of accuracy in the 
solution (this is discussed in subsequent paragraphs of the section), and 
the drag factor matrix [~ K DL ~J is made to be proportional to the mass matrix 
['Mj, then Eq. 3.40 can be written as 

FMj{x $ } + [C s ]{x $ } + ar-Mj{x s > + [K]{x $ } = [K]{x g > (3.41) 

In order to transform Eq. 3.41 into a "normal" coordinate system, the trans- 
formation 



(x s > = [<f>]{Q} (3.42) 

is made. Substituting 3.42 into Eq. 3.41 and premultiplying the result by 

[ 4 >] T , 

D*M*J{Q} + r'C*J{Q} + arM*J{Q} + FK*J{Q} 

= [4>] T [K]{x g } (3.43) 

Define 

{y} = [<f>] T [K] (I } (3.44) 

and Eq. 3.43 can be written as 

+ rC*J{Q} + + C'K*J{Q} = { y}x n 

^ y 



(3.45) 
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Equation 3.45 represents nm uncoupled single degree of freedom modal equations 
of motion. The equation of motion for the ith mode is 



M* . Q . + C* Q. + <xM*.Q. + K*,Q. = y.x 

n i s ii l li i n i l g 



(3.46) 



or, 

«i + 26 s “i«i + 26 h “11«1 + “ii Q i = <tV ><“f i x g> < 3 - 47 > 
i i n ii w ii 

Since the left side of Eq. 3.47 is linear with respect to the unknown Q, the 
response quantities of interest in the ith mode can be obtained by computing 
the desired response quantities from 



9 9 

q + 26 s wq + 2(3pU>q + o> q = to Xg (3.48) 

and multiplying the result times the ith mode participation factor for hydro- 
dynamic parameters . The participation factor for hydrodynamic parameters 
can be seen from Eqs. 3.44 and 3.47 to be 



{ HDPF} = r'K*-J" 1 [<D] T [K]{I} 



(3.49) 



Once the modal quantities are obtained in the Q coordinate system, Eq. 3.42 
can be used to obtain an approximation of the actual hydrodynamic parameter 
of interest at each node in much the same manner as an estimate of the maxi- 
mum deformation of each node is obtained by modal superposition. 

To answer the first question raised in this section, the suitability 
of dropping the velocity force, [C s ]{x g }, from Eq. 3.40, response calculations 
were made for the following two equations of motion: 

C'Mj{x s > + [C s ]{x s > + rK DL -J{x s } + [K]{x s } 

= [K](x g } + [C s ]{x g > 



(3.40) 



75 



and 

['Mj{x s > + [C s ]{x s > + rK DL -J{x s } + [K]{x s > = [K]{x g } (3.50) 

The quantities ^ l x s l max } > ^ < l x s l > >> { l^sLax >5 and for systems 1 ’ 

4, and 7 in the mid-frequency range were calculated for the El Centro earth- 
quake from both Eqs. 3.40 and 3.50. The damping matrices [C g ] used in the 
series of calculations were such that they yielded damping ratios of .02, 

.05, and .10 in each mode and IT~ waS set t0 zero * For each structural 
system and value of [C ], the hydrodynamic parameters computed from Eq. 3.40 
were within 5 percent of the corresponding response quantities computed from 
Eq. 3.50. The results of the calculations for system 4 and a damping matrix 
which yields .02 in each mode are shown in Table 3.34. In another set of 
calculations, the damping matrices [C $ ] were made to be proportional to the 

mass matrix and the numerical values of [C ] were such that the first mode 

s 

damping ratios were .02, .05, and .10. Again, for each structural system 

and value of [C s ], the response quantities computed from Eq. 3.40 were within 

5 percent of the corresponding response quantities computed from Eq. 3.50. 

The numerical results were such that no clear trend resulted that indicated 

which equation yielded consistently higher or lower results. The hydrodynamic 

parameters were computed for system 4M in the mid-frequency range, with a 

structural damping matrix which yielded modal damping ratios of .02 in each 

mode, for the El Centro earthquake. The drag factor matrix was linearized by 

using the ground velocity characteristics and made to be proportional to the 

mass matrix. In one set of computations, the velocity force [C ]{x } in Eq. 

^ y 

3.40 was included and in another set of computations, the velocity force 
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[C ]{x } was neglected. The two sets of numerical results differed by less 
s g 

than 5 percent. 

From these calculations, it is concluded that the velocity force 
[C $ ] {Xg} can be dropped from Eq. 3.40 without sacrificing the accuracy of 
the results by more than about 5 percent, for values of [C s ] which yield 
modal damping ratios of up to about .10. 

In order to show the similarity of form between Eq. 3.48 and the 
linear equation of motion for the single degree of freedom oscillator dis- 
cussed in Chapter 2, it is instructive to recall Eq. 2.7, repeated here for 
convenience. 



Z + 26„wZ + 28t,wZ + a) Z = -x - 28.wx 
s h g h c 



(3.51) 



Since Z = x - x , Eq. 3.51 can be written 
s g 

9 9 

x„ + 26 u)X o + 26 ^ 0 )^ + co x„ = co x„ + 26 cox 

s s s hs s g sg 



(3.52) 



It was seen in Chapter 2 that, for values of 3 tota -j up to about .10, the 
hydrodynamic parameters computed from Eq. 3.51 were dependent upon the amount 
of total damping in the system (6 $ + 6^) and not on the distribution of the 
damping between structural and hydraulic. If all of the damping in Eq. 3.51 
is considered to be structural and all of the damping in Eq. 3.52 is considered 
to be hydraulic, the velocity force on the right side of both equations can 
be dropped, and Eqs. 3.51 and 3.52 become 



and 



Z + + in Z = -x 

t g 



x s + 26 t u)X s + 0 ) 2 x s = w 2 x g 



(3.51a) 



(3.52a) 
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In other words, for values of 3 t up to about .10, the hydrodynamic parameters 
lx I ,<|x|>, Ixl, and <|x |> computed from Eq. 3.51a are within 2 
to 4 percent of those which would be computed from Eq. 3.52a for the same 
system and base disturbance. 

It is seen that Eq. 3.48 is of the same form as Eq. 3.52a, so if 
the hydrodynamic parameters can be calculated, or estimated (as from an ideal- 
ized response spectrum) for Eq. 3.52a with satisfactory accuracy, then modal 
hydrodynamic parameters |q| x> < l ( ll > > I ^ I max ’ anc * < ^l > can a ^ so * 3e ca ^ cu l at ed 
or estimated. Once the modal response quantities I q I max » < hl > ’ |q| max » anc * 

<|q|> are obtained and multiplied times the appropriate modal hydrodynamic 

• • 

parameter participation factor to obtain |Ql max » <|Q|>> I ^ I max ’ anc * < I Q I > > 

the problem arises of how to combine the modal responses. 

To obtain an estimate of {lx | „ } and {|x | „ }, the modal responses 

1 s' max 1 s' max 

can be combined in any one of several methods as is the case for obtaining 
{|Z| m ax^ An upper bound is obtained by summing the absolute values of the 
modal responses (ZABS). 

An upper bound on the quantities { < I x s I >> anc * {< | x $ | >} can also be 
obtained by the ZABS method of modal combination, but this is not readily 
apparent. That the ZABS method of modal combination yields an upper bound 
on the quantity {< | x s | >} (or { < | x s | >} ) is shown as follows: 

{x s > = [<f>]{QJ (3.42) 

Assume that the vector (Q ( t ) } has been obtained, then, at the jth node of 
the system, 

x (t) - ♦j 1 Q 1 (t) + ■tj 2 q 2 (t) + • • • ♦j.jwo 

J 



(3.53) 
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Taking the absolute value of both sides of Eq. 3.53 

I x s - (t) I = l^jlV 1 ^ + + ' • • ^j.nmW 1 ^ ^ 3 ' 54 ^ 

J 

An upper bound on |x (t)| is seen from Eq. 3.54 to be 

j 

|x $ (t)| < Uj-j I lQ-| (t) I + I4»j2 1 l Q 2 (t) l + * • * l Q nm (t) l 

(3.55) 

Integrating both sides of Eq. 3.55 over the time of response. 



1 

t 



t 

)x ( t ) | d t 

o j 




r t 

0 



1 4*j2 1 1^2^^ I 



+ . 




IIQ nra (t)|]dt 



(3.56) 



or, 



<|xj> 5 l4*j-| l<lQ-i l> + 

J 



+ . . . U. I < I Q |> 
IY j,nm' |V nm' 



(3.57) 



Which is seen to be the ZABS method of modal combination. 

Rules for constructing idealized hydrodynamic parameter response 
spectra from which |q| , <|q|>, |q| , and <|q|> may be estimated for each 

maX lllaX 

• • 

mode (and then |Q| m , < | Q | > * lQl max » anc * < |Q| > are computed from the q's and 
the hydrodynamic parameter participation factors) of the multi-degree of 
freedom system are discussed in the next section. 
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3.5 Idealized Elastic Hydrodynamic Parameter Response Spectra 



Idealized elastic hydrodynamic parameter response spectra can be 
constructed for each mode of the multi-degree of freedom system in the fol- 
lowing manner: 

1. Compute the natural frequencies, ['wj, and mode shapes, [<J>], 
for the system. 

2. Establish the structural damping ratio, 6 , for each mode. 

s i 

3. Compute and HDPF.. from 



a = 



/ nn 

I k d.. 
3=1 .13 

nn 

1 M. . 
3=1 JJ 



\ 


• 2 r i \ 
<X 9 |X 9 l> 




•2 




<x > 


/ 


9 / 



(3.22) 



r2$ h wj = a["I-J 

{HDPF} = r K*J" 1 [^] T [K]fI> 



(3.24) 

(3.49) 



4. With the quantities u>.., 3 S . , B h ., x g max (or <|x g |>), x g max 
(or <|x |>), and x (or <|x |>), construct the idealized 

elastic hydrodynamic parameter response spectrum for the ith 
mode exactly as described in Section 2.7 for the single degree 
of freedom system. 

When an idealized spectrum has been drawn for each mode of the sys- 
tem, maximum modal hydrodynamic parameters are obtained from the spectra, 
multiplied times the appropriate modal hydrodynamic parameter participation 
factor , and combined in some manner to yield an approximate solution of the 
hydrodynamic parameter of interest. 
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CHAPTER 4 
APPLICATIONS 



4.1 General 

In Chapters 2 and 3 a method was developed for linearizing the non- 
linear hydrodynamic drag force on a single degree of freedom oscillator and 
a multi-degree of freedom lumped parameter model. The method has been shown 
to yield satisfactory results for systems subjected to earthquake excitation 
of the type used in this thesis. Although the method of linearizing employed 
in the multi-degree of freedom system places a linear damper at the level of 
the top mass where a nonlinear damper does not exist, it has been shown to 
yield maximum response quantities which are within a very few percent of the 
response quantities computed for the nonlinear system. Further, linearizing 
the drag force by using velocity characteristics of the ground motion has been 
shown to yield generally conservative results in all but the low frequency 
range. The linearization technique has been developed in such a manner that 
the resulting linear simultaneous differential equations of motion for the 
multi-degree of freedom system possess classical normal modes of response. 

The existence of classical normal modes permits the use of normal mode 
superposition, with the modal response quantities determined by response 
spectrum techniques, as an approximate method of computing the response of 
multi -degree of freedom systems. The approximate method can be used for dis- 
placement response quantities and hydrodynamic parameter response quantities. 
The research reported in this thesis is not intended to reinvent the method of 
solution by the superposition of normal mode solutions as determined by 
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response spectra, but is devoted primarily to developing the linearizing 
technique which insures the existence of classical normal modes and yields 
generally conservative results. 

This chapter demonstrates the use of the principles and methods 
developed and discussed in Chapters 2 and 3 in computing the response of 
a series of actual offshore tower designs and one of the test systems des- 
cribed in Chapter 3. 

4.2 Displacement Response 

Table 4.1 lists the properties of a series of offshore towers de- 
signed for water depths of 400, 600, 800, and 1000 feet. 4 A structural damp- 
ing matrix [C s ] which yields modal damping ratios of .01 in each mode is com- 
puted for each tower as described in Chapter 3. The structures are subjected 
to the El Centro earthquake depicted in Fig. 1.2, and various response quan- 
tities are computed. In the first set of computations, the equations of mo- 
tion are numerically integrated and response quantities are found for the 
case in which there are nno hydrodynamic drag forces acting on the structure. 

The same response quantities are next computed for the case in which the actual 
nonlinear drag forces are acting on the structure. Finally, the nonlinear 
drag forces are linearized by using the velocity characteristics of the 
ground motion, and the response quantities are computed by the combination 
(EABS and SRSS) of maximum modal responses found from the true response spec- 
trum for each of the first five modes of the structures. The amount of linear 
hydraulic damping in each mode found from the linearization procedure is shown 
in Table 4.2. The response quantities of interest in these calculations are 
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the maximum displacement between the top deck and the ground, the maximum over- 
turning moment, and the maximum base shear expressed as a percent of the total 
tower weight, including the added mass of the structure. The results of the 
calculations are shown in Fig. 4.1. 

A study of Fig. 4.1 reveals that, for this series of towers sub- 
jected to the El Centro earthquake, the nonlinear hydrodynamic drag forces 
reduce the response of the tallest tower in the top deck displacement and 
overturning moment quantities by approximately 30 percent from the "no drag" 
case. The maximum base shear for the tallest tower is reduced only about 8 
percent by the nonlinear hydrodynamic drag forces. Exact response calculations 
were made for the other three earthquake records used in this thesis and the 
reduction in response due to drag was found to be somewhat lower in magnitude 
than for the El Centro earthquake. For the tower in 1000 feet of water, the 
reduction (due to hydrodynamic drag) in the maximum top deck displacement 
relative to the ground was found to be about 8.5 percent for the Taft earth- 
quake, 11.5 percent for the Ferndale earthquake, and 17 percent for the San 
Fernando earthquake. In general, the taller, more flexible towers showed 
more reduction in response due to the hydraulic drag than did the shorter, 
stiffer towers. No generalizations can be made, however, regarding the mag- 
nitude of the reduction in response relative to tower height or stiffness 
since the actual response is such a complex combination of variables in both 
the structure and the earthquake. 

In the superposition of maximum modal responses computed from the 
true response spectrum for each of the first five modes of vibration, it was 
found that the higher modes contributed rather strongly to the maximum base 
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shear. The contribution of each mode to the maximum base shear computed by 
the EABS combination of modal responses is shown in Table 4.3. Although the 
SRSS combination of modal responses appears from Fig. 4.1 to provide a very 
good estimate of the response of this series of towers to the El Centro earth- 
quake, generalities concerning the accuracy of one method of modal combination 
relative to another method of modal combination cannot be made. 13 

The response of the towers is computed by the combination of maxi- 
mum modal responses (ZABS and SRSS) found from the idealized deformation re- 
sponse spectrum for each of the first five modes of vibration. The idealized 
spectra were constructed in accordance with the procedure described in Chapter 
3 for ground motion quantities equal to those of the El Centro earthquake. The 
idealized spectra for the first five modes of tower C are shown in Fig. 4.2. 

The estimated hydraulic damping ratios in the 3rd, 4th, and 5th modes of the 
structure are such that the resulting idealized spectra for those modes cannot 
be shown separately and distinctly on the scale of Fig. 4.2. It is noted that, 
for the idealized spectra shown in Fig. 4.2, the frequencies of the system are 
well removed from the high frequency tail region of their respective response 
spectrum. This was found to be the case for each of the towers A, B, C, and 
D. The results of the computations using the idealized spectra are shown in 
Fig. 4.3. 

4.3 Hydrodynamic Parameters 

To demonstrate the computation of the hydrodynamic parameters 

|x s l max , <|x s l>. I^ s lmax’ anc * < l x s^ > at eac ^ noc * e a mu lti-degree °f freedom 
system, the five degree of freedom system 4M is used. The properties of this 
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system, with a structural damping matrix that yields a damping ratio of .02 
in each mode, are shown in Table 4.4. The magnitude of [K] was chosen to 
yield frequencies of the system that fall into the velocity amplification 
region of the spectrum. The hydrodynamic parameter response quantities are 
computed by the numerical integration of the equations of motion and by the 
combination (ZABS and SRSS) of maximum modal responses found from the true 
spectral response in each mode. The results of the computations are shown 
in Table 4.5. It is seen from Table 4.5 that the combination of modal re- 
sponses provides results which are, in general, sufficiently close to the 
exact response at each node to enable the selection of inertia and drag 
coefficients Cj and C Q with reasonable confidence. In lieu of using the 
true spectral response quantities, which require the use of a digital com- 
puter to obtain for each mode, an idealized spectrum can be constructed for 
each mode in the manner set forth in Chapter 3. 
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CHAPTER 5 

CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 



5.1 Conclusions 



In this thesis, an approximate method has been developed for com- 
puting the response of a multi -degree of freedom system, partially submerged 
in water, to earthquake excitation. In Chapter 2, the effects of the non- 
linear hydrodynamic drag force on a single degree of freedom oscillator are 
shown, and a generally conservative method of linearizing the drag force is 
developed. Amplification factors are developed which permit the construction 
of idealized response spectra for deformation and hydrodynamic parameters. 
Chapter 3 extends the principles from the single degree of freedom oscillator 
to multi-degree of freedom systems and develops a generally conservative 
method of linearizing the drag forces. The method of linearizing insures the 
existence of classical normal modes in the equations of motion, which in turn 
permits the response quantities of interest to be obtained by the superposition 
of normal mode responses. In Chapter 4 the response to the El Centro earthquake 
of a series of offshore towers is computed using the methods developed in 
Chapter 3. 



Several conclusions may be drawn from the development of principles 

and methods and the results of sample problems in the preceding three chapters. 

1. The maximum relative displacement response of a single degree 

of freedom oscillator, with drag characteristics described by 
K D *a max 

the parameter ^ , approaches a constant value at high 

frequencies as a result of the nonlinear hydrodynamic drag 
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force growing very large relative to the inertia force in the 
equation of motion. The magnitude of this relatively constant 
high frequency tail can be estimated for earthquake motions by 
the expression 



K n x 1 


/ X 2 


D g max 


g max 


K 


X 


1 


g max 



2. Linearizing the drag force on a single degree of freedom oscil- 
lator by using the absolute velocity response of the oscillator 
yields slightly conservative results for the maximum response 
to earthquake motions until the high frequency, high drag tail 
begins. 

3. Linearizing the drag force on a single degree of freedom oscil- 
lator by using the velocity characteristics of the ground motion 
yields conservative results for the maximum response to earth- 
quake motions except for low frequency systems and very high 
frequency systems. 

4. The conservatism introduced by using the ground motion velocity 
characteristics in the linearization of the drag force is re- 
duced as the amount of structural damping in the system is in- 
creased. 

5. The maximum response of the single degree of freedom oscillator 
to earthquake motions is sensitive to the total amount of linear, 
velocity dependent damping in the system, and not to the distri- 
bution of this damping between structural and hydraulic, up to 
linear hydraulic damping ratios of 10 percent of critical. 
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6. The nonlinear hydrodynamic drag has very little effect on the 
elasto-plastic deformation response spectrum in the low fre- 
quency region. 

7. A high frequency tail exists in the elasto-plastic deformation 
response spectrum for a system with nonlinear hydrodynamic drag 
forces. Large amounts of drag cause greater reductions in the 
yield deformation, Z^, from the maximum elastic deformation for 
a given ductility factor than do small amounts of drag. 

8. The nonlinear hydrodynamic drag force has a damping effect on 
the elasto-plastic deformation response spectrum in the mid- 
frequency range. The damping effect of the drag becomes less 
pronounced as the ductility factor increases. 

9. Idealized response spectra which consider the effects of fluid- 
structure interaction can be constructed as described in this 
thesis. 

10. The linearization of the drag forces on a multi-degree of freedom 

system by using the absolute velocity response at each node and 
forcing the linear drag factor matrix to be proportional 

to the mass matrix t'Mj as described herein is satisfactory. 

11. Linearizing the drag forces on a multi-degree of freedom system 
by using the velocity characteristics of the ground motion and 
forcing the drag factor matrix to be proportional to the mass 
matrix as described herein yields generally conservative results 
except for low frequency systems. 

12. For the multi -degree of freedom systems studied in this thesis. 
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the response to earthquake excitation is relatively insensitive 
to the spatial distribution of the linear drag forces. 

13. The hydrodynamic drag forces can cause a significant reduction 
in the response of lightly damped systems to earthquake motions. 

14. The amount of reduction in response due to hydrodynamic drag 
varies considerably from one response quantity to another and 
from one earthquake to another. 

15. The reduction in response due to hydrodynamic drag is decreased 
as the amount of structural damping in the system is increased. 

16. Higher modes of a multi-degree of freedom system can contribute 
appreciably to the total response of the system. 

5.2 Recommendations for Further Study 

The approximate methods of analysis developed in this thesis are 
limited by the types of multi -degree of freedom systems used in the investi- 
gation and by the inherent limitations in the hydrodynamic force equation 
used. Further, the study in this thesis is devoted to excitation of the 
structure in still water, and the ocean is never completely at rest. There 
are several areas in which further studies should be made of the response of 
offshore structures to earthquake excitation. 

1. The validity of the lumped parameter model should be tested 
through the use of field measurements of the response of actual 
offshore structures to strong motion earthquakes. 

2. Model tests of single degree of freedom oscillators and multi - 
degree of freedom structures should be conducted for the types 
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of base motion considered in this thesis. Such model tests will 
be useful in verifying the principles developed herein. 

3. Since the ocean is never actually at rest, analytical and ex- 
perimental studies should be conducted on the response of sub- 
merged and partially submerged structures subjected to earthquake 
motions combined with forces resulting from waves and currents. 

4. The studies in this thesis assume ground motion in one direction 
only. It is recommended that studies be made of the response 
when the two orthogonal components of the horizontal ground motion 
are acting on the structure at the same time. 

5. The earthquake response of the structure in torsion and in the 
vertical direction should be investigated. 

6. It is recommended that the effects of foundation and soil con- 
ditions be included in the model. There are several reasonable 
approaches to modeling the foundation which should be investi- 
gated. 23-25 

7. The hydrodynamic force equation used in this thesis (Mori son's 
equation) is based on studies with relatively rigid piles. There 
is some evidence 19 * 20 that the forces on flexible, oscillating 
structures may be better expressed through the use of different 
drag and inertia coefficients, or even through the use of a 
different force expression. Much more information is needed 
concerning the hydrodynamic forces acting on flexible, oscillating 
structures. 

8. The linearizing technique should be tested for suitability of 
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use when the structure has radically different distribution of 
mass, stiffness, and drag properties from those test structures 
used in this thesis. 

9. It is recommended that the effects of the hydrodynamic forces 
on the inelastic response of the structure be studied in greater 
depth. While design for elastic response is satisfactory for 
earthquakes which may be expected to occur frequently, elastic 
design for infrequent, severe earthquakes may result in ex- 
cessively conservative designs. Some response into the in- 
elastic range should be permitted for the occurrence of an un- 
usually strong earthquake, and approximate methods are desirable 
for the estimation of this inelastic response. 

10. It is recommended that approximate methods be developed for 
the estimation of the response to stochastic representations 
of the forces resulting from earthquakes, waves, and currents. 
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Table 2.1 

Characteristics of Series Representation 
of Base Acceleration 



Quantity Maximum Amplitude Time of Occurrence 



x 2.56 g 0.35 1 1 

g d 

x 3.85 ft/sec 0.42 t, 

9 a 

x 0.20 ft 0.50 t. 

g a 

Xg = -0.20401 cos 2frt - 6.87108 cos 4irt + 38.35016 cos 6iTt 
- 41.35223 cos 8irt - 2.84232 cos 10rrt 
+ 12.91937 cos 1 2tt t 

(Xg in ft/sec^) 



Comparison of Estimated | Z | max With |Z| max Computed at 
High Frequencies by Numerical Integration 
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Table 2.3 

Distribution of Total Damping Between Structural 
and Hydraulic for Comparison Purposes 



^total 




3 h For 






B $ = 0 


3 $ = .02 


3 $ = .05 


3 S = .10 


.02 


.02 


0 






.05 


.05 




0 




.07 




.05 


.02 




.10 


.10 




.05 


0 


.12 




.10 




.02 


.15 






.10 


.05 



Table 2.4 

Elastic Deformation Spectrum Amplification Factors 



W*> 


Ampl ification 


Factor 


Displacement 


Velocity 


Accel eration 


0 


2.5 


4.0 


6.4 


0.5 


2.2 


3.6 


5.8 


1 


2.0 


3.2 


5.2 


2 


1.8 


2.8 


4.3 


5 


1.4 


1.9 


2.6 


7 


1.2 


1.5 


1.9 


10 


1.1 


1.3 


1.5 


20 


1.0 


1.1 


1.2 
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Table 2.5 

Elastic Hydrodynamic Parameter Spectrum 
Amplification Factors 



l x sLax Spectrum 
Amplification Factors 



total 

(*) 


Displacement 


Velocity 


0 


2.4 


5.2 


1 


2.2 


4.4 


2 


2.1 


3.6 


5 


1.8 


3.1 


10 


1.6 


2.1 




<|x |> Spectrum 




Ampl ification 


Factors 


'total 

(%) 


Displacement 


Velocity 


0 


3.8 


10.9 


1 


3.4 


8.5 


2 


3.0 


6.8 


5 


2.3 


4.2 


10 


1.8 


3.3 
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Table 2, 


.5 (Continued) 




0 total 

(%) 


l*sLx S » ect ™" 
Amplification Factors 


Displacement 


Velocity 


Acceleration 


0 


2.2 


4.1 


7.1 


1 


1.9 


3.1 


4.9 


2 


1.6 


2.6 


3.7 


5 


1.4 


2.0 


3.1 


10 


1.2 


1.8 


2.6 






< | x | > Spectrum 




^total 

(%) 




Amplification Factors 


Displacement 


Velocity 


Acceleration 


0 


3.4 


10.1 


13.5 


1 


3.0 


6.3 


6.2 


2 


2.7 


5.1 


4.3 


5 


2.1 


3.0 


2.9 


10 


1.6 


2.2 


2.2 
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Table 2.6 



Actual and Approximate Reduction Factors for the 
High Frequency Region of Inelastic Response 
Spectra. El Centro Earthquake 

. 002> K D X 9 max . 050 
S l\ 



Ductility Factor 

y 


Z 

y 

(inches) 


i z U 

(inches) 


Actual 

Reduction Factor 
for Z y 


y 


2y - 1 


1 


0.63676 


0.63676 


1.00 


1.00 


1.5 


0.49526 


0.74289 


0.77 


0.75 


2 


0.42835 


0.85670 


0.67 


0.67 


3 


0.35975 


1.07925 


0.56 


0.60 


4 


0.31523 


1.26093 


0.49 


0.57 



Note: All calculations are at a frequency of 10 cps. 



Relative Values of Mass, Stiffness, and Drag Factor for 5 Degree of Freedom Systems 
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CASE 2 - GROUND VELOCITY IS LSEO TO LINEARIZE THE DRAG FORCE 



Response of System 1L to El Centro Earthquake, 1940, N-S Component. 
Direct Integration Solution. 3 = .05 
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Response of System 1M to El Centro Earthquake, 1940, N-S Component. 
Direct Integration Solution, g = .01 
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CASE 2 - GROUND VELOCITY IS USEC TO LINEARIZE THE DRAG FORCE 



Response of System 2M to El Centro Earthquake, 1940, N-S Component. 
Direct Integration Solution. 3 = .05 
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CASE 2 - GROUND VELOCITY IS USED TO LINEARIZE THE DRAG FORCE 



Response of System 3L to El Centro Earthquake, 1940, N-S Component. 
Direct Integration Solution. 3 = -01 
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CASE 1 - VELOCITY OF EACH MASS (COMPUTED IN ACTUAL DRAG PR08LEM) IS USED TO LINEARIZE THE DRAG FORCE 
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CASE 2 - GRCUND VELOCITY IS LSEC TO LINEARIZE THE DRAG FORCE 



Response of System 6M to El Centro Earthquake, 1940, N-S Component. 
Direct Integration Solution. 6 = .01 
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Response of System 7M to El Centro Earthquake, 1940, N-S Component. 
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Response of System 7M to El Centro Earthquake, 1940, N-S Component. 
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Table 3.34 



Hydrodynamic Parameters for System 4, Mid-Frequency Range. 
El Centro Earthquake. $ $ = .02 in Each Mode 



Parameter 


Node 


Including 

CC s ]{x g } 


Excluding 

[ y<y 


^ Imax^ 


1 


18.95250 


18.98743 


(inches) 


2 


16.31281 


16.37181 




3 


13.45125 


13.47894 




4 


12.30791 


12.21299 




5 


10.53108 


10.50990 


{<|x s |>} 


1 


6.28175 


6.26993 


(inches) 


2 


5.21131 


5.20904 




3 


4.48454 


4.48830 




4 


3.75600 


3.75798 




5 


2.81634 


2.81841 


^ l x s Imax ^ 


1 


32.53683 


32.50833 


(in. /sec) 


2 


23.84686 


24.27509 




3 


19.37517 


19.68547 




4 


21.22523 


21.28960 




5 


15.38559 


15.47496 


«|x s |» 


1 


10.26027 


10.26663 


(in. /sec) 


2 


7.95916 


7.94259 
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6.83778 


6.82866 
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5.61321 


5.60034 
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Table 4.2 

Estimated Linear Hydraulic Damping Ratios 
for Towers A, B, C, and D 



Mode 




Towers 




A 


B 


C 


D 


1 


0.00587 


0.00632 


0.01178 


0.01015 


2 


0.00251 


0.00304 


0.00576 


0.00541 


3 


0.00145 


0.00168 


0.00315 


0.00328 


4 


0.00107 


0.00126 


0.00246 


0.00233 


5 


0.00085 


0.00115 


0.00199 


0.00183 



Table 4.3 

Contributions of Each Mode to the Maximum 
Base Shear. ZABS Modal Combination 




Contribution (as 


% of Total) of Mode 




Tower 


1 2 


3 


4 


5 


A 


32.8 40.2 


18.7 


6.8 


1.5 


B 


33.0 28.3 


24.0 


13.6 


1.1 


C 


19.3 24.6 


15.6 


35.0 


5.5 


D 


38.9 41.8 


17.9 


1.4 


0.0 
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Table 4.4 

Properties of System 4M, Mid-Frequency Range. 
3 s = .02 in Each Mode 



[~MJ: 


1 1 


1 


2 


3 


2 

lb-sec /in. 


rK 0 J: 


0 .00457 


.00457 


.00914 


.01371 


lb-sec /in. 


[K]: 


15.32 -15.32 


0 


0 


0 






-15.32 45.96 


-30.64 


0 


0 








61.29 


-30.64 


0 


lb/in. 




Sym. 




61.29 


-30.64 












76.61 




[C s ]: 


0.13917 -0.06907 


-0.01734 


-0.01149 


-0.00479 






0.24264 


-0.09762 


-0.02638 


-0.00830 








0.28693 


-0.10783 


-0.01633 


Ib-sec/in. 




Sym. 




0.40007 


-0.12933 












0.58443 




Modal Properties 
















Hydrodynamic Parameter 


Mode 


Frequency 3^ 


(estimated ) Parti ci pati on 


Factor 




(cps) 


(none) 




(none) 




1 


0.24705 


0.00794 




2.49115 




2 


0.55706 


0.00352 




1.13695 




3 


0.84337 


0.00233 




0.66289 




4 


1.06231 


0.00185 




0.25179 




5 


1.52207 


0.00129 




0.01573 
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Table 4.4 (Continued) 



Spectral Response Quantities (True Spectra) 



Mode 


l q Lax 


<hl> 


|q | 

Ml max 


A 

• cr 

V 




(in. ) 


(in. ) 


(in. /sec) 


(in. /sec) 


1 


12.618 


4.044 


16.648 


5.295 


2 


13.382 


3.197 


22.807 


7.056 


3 


9.316 


3.096 


25.950 


10.182 


4 


10.376 


2.926 


38.248 


9.798 


5 


10.811 


2.712 


31.994 


9.387 
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Table 4.5 

Hydrodynamic Parameters for System 4M. 
El Centro Earthquake 



Response 

Quantity 


Node 


Exact 

Value 


Modal Superposition - 
Using True Spectral Response 


EABS 


SRSS 




1 


17.82343 


31 .32392 


21.00747 




2 


16.05878 


21.29961 


16.04786 


{|x sLx } 


3 


13.26795 


18.78246 


13.53510 


(inches) 


4 


11.71001 


17.13259 


11.14007 




5 


10.23945 


11.85761 


6.68999 




1 


5.09378 


9.28465 


6.43887 




2 


4.08737 


6.63971 


5.12691 


«|x s l» 


3 


3.44525 


5.81354 


4.30800 


(inches) 


4 


2.63410 


4.96109 


3.32816 




5 


1.45048 


3.42885 


1.91468 




1 


31.62296 


50.63429 


30.31177 




2 


22.94491 


35.58812 


22.48838 


‘l^sLax 1 


3 


17.93618 


30.87149 


19.23872 


(in. /sec) 


4 


19.32336 


26.62459 


16.66792 




5 


13.95569 


22.63737 


12.06787 




1 


8.98920 


16.32071 


9.69329 




2 


6.82424 


11.41066 


7.21779 


{<|x s l>) 


3 


5.89449 


10.20436 


6.32494 


(in. /sec) 


4 


4.96801 


8.89230 


5.18860 




5 


4.10293 


7.49280 


4.08744 



Operating 

Platform 
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100.0 
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33S/NI ‘UI3013A 



Fig. 1.2 Elastic Deformation Response Spectra for El Centro Earthquake 



DISPLACEMENT (INCHES) VELOCITY (INCHES/SECOND) PCCE-ERAT I ON (C) 

-10.0 -S.O 0.0 S.O -7.0 0.0 7.0 -0.20 0.0 0.20 
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Fig. 1.3 Ground Acceleration, Velocity, and Displacement of 
El Centro Earthquake, 1940, N-S Component 



DISPLACEMENT (INCHES) VELOCITY (INCHES/SECOND) ACCELERATION (G) 

- 2.0 0.0 2.0 - 4.0 0.0 4.0 - 0. 10 0.0 0.10 
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Fig. 1.4 Ground Acceleration, Velocity, and Displacement of 
Taft Earthquake, 1952, N-21-E Component 



DISPLACEMENT (INCHES) VELOCITY (INCHES/SECOND) ACCELERATION IG) 

- 5,0 0.0 5.0 - 9.0 0.0 9.0 - 0.10 0.0 0.10 
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0.0 3.0 6.0 9.0 12.0 15.0 18.0 21.0 2M.0 27.0 30.0 



TIME (SECONDS) 

Fig. 1.5 Ground Acceleration, Velocity, and Displacement of 
Ferndale Earthquake, 1954, N-44-E Component 



DISPLACEMENT (INCHES) VELOCITY (INCHES/SECOND) ACCELERATION (G) 

-15.0 0.0 15.0 -2S.0 0.0 25. 0 -1.0 0.0 
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0.0 2.0 1.0 6.0 8.0 10.0 12.0 1M.0 16.0 18.0 20.0 



TIME (SECONDS) 

Fig. 1.6 Ground Acceleration, Velocity, and Displacement of San 
Fernando Earthquake (Pacoima Dam Record), 1971, 

S-16-E Component 



0001 
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Fig. 2.1 Elastic Deformation Response Spectra for x = 16.1 sin 2 ir t (ft/sec^). 6 = 0.01 
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Fig. 2.2 Elastic Deformation Response Spectra for x_ = 16.1 sin 2 ir t (ft/sec 2 ). 6 = 0.05 
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Fig. 2.3 Elastic Deformation Response Spectra for Six Term Series Base Acceleration Record. 3 =0.01 
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Fig. 2.5 Elastic Deformation Response Spectra. 
El Centro Earthquake. 6 = .00 
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Elastic Deformation Response Spectra 
El Centro Earthquake, B = .02 
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Elastic Deformation Response Spectra. 
El Centro Earthquake. B = .05 
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Fig. 2.8 Comparison of Elastic Deformation Response Spectra With Linear and Nonlinear Drag 

for Series Representation of Base Disturbance. 8 = 0.01, KJ< /K = 0.20 
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Fig. 2.9 Comparison of Elastic Deformation Response Spectra With Linear and Nonlinear Drag 
for Series Representation of Base Disturbance. 6 = 0.05. ICx /K = 0 20 
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Fig. 2.10 Comparison of Elastic Deformation Response Spectra With Linear and Nonlinear Drag 
for El Centro Earthquake. 8 S = 0.02, max /K = 0.50 
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Fig. 2.11 Comparison of Linear and Nonlinear Solution. (Linearized 
by Oscillator Velocity) El Centro Earthquake. 6 = .00 
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Fig. 2.12 Comparison of Linear and Nonlinear Solution. (Linearized by 
Oscillator Velocity) El Centro Earthquake. 6 = .02 
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Fig. 2.13 Comparison of Linear and Nonlinear Solution. (Linearized 
by Oscillator Velocity) El Centro Earthquake. 6 = .05 
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Fig. 2.14 Linearizing Parameter for Single Degree of Freedom 
Oscillator. El Centro Earthquake. 3 = .00 
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Fig. 2.15 Linearizing Parameter for Single Degree of Freedom 
Oscillator. El Centro Earthquake. 3 = .02 
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Fig. 2.16 Linearizing Parameter for Single Degree of Freedom 
Oscillator. El Centro Earthquake. 6 = .05 
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Fig. 2.17 Comparison of Nonlinear Drag Spectrum With Linear Drag Spectrum (Linearized 

by Ground Velocity). El Centro Earthquake. 3 = 0.02, K n x '/K = 0.006 
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Fig. 2.18 Comparison of Nonlinear Drag Spectrum With Linear Drag Spectrum (Linearized 
by Ground Velocity). El Centro Earthquake, g_ = 0.02, K n x„ mav/ /K = 0.08 
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by Ground Velocity). El Centro Earthquake. 
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Fig. 2.20 Comparison of Linear and Nonlinear Solution (Linearized 
by Ground Velocity). El Centro Earthquake. 6 =0.02 
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Fig. 2.21 Elastic Deformation Response Spectra for El Centro Earthquake 
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Fig. 2.21 (Continued) 
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Fig. 2.22 Elastic Deformation Response Spectra for Taft Earthquake 
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Fig. 2.22 (Continued) 
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Fig. 2.23 Elastic Deformation Response Spectra for Ferndale Earthquake 
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(d) B $ - 10% 6 h = 0, 1, 2, 5, 10% 



Fig. 2.23 (Continued) 
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Fig. 2.24 Elastic Deformation Response Spectra for San Fernando Earthquake 
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(d) B s = 10* B h = 0, 1, 2, 5, 10* 



Fig. 2.24 (Continued) 
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Fig. 2.25 Comparison of Spectra for Two Distributions of Total 
Damping Ratio. Ferndale Earthquake. 3 tota -j = .07 
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Fig. 2.26 Construction of Idealized Elastic Deformation Response Spectrum 
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Fig. 2.27 Elastic Hydrodynamic Parameter Response Spectra. 

El Centro Earthquake. 3 total = 0, 1, 2, 5, 10% 
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Fig. 2.27 (Continued) 
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Elastic Hydrodynamic Parameter Response Spectra. 
Taft Earthquake. 3 total = 0, 1, 2, 5, 10% 
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Fig. 2.28 (Continued) 
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Fig. 2.29 Elastic Hydrodynamic Parameter Response Spectra. 

Ferndale Earthquake. £ tota -] = 0> 1» 2, 5, 10% 
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Fig. 2.29 (Continued) 
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Fig. 2.30 Elastic Hydrodynamic Parameter Response Spectra. 

San Fernando Earthquake. & total = 0, 1, 2, 5, 10% 
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Fig. 2.30 (Continued) 
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Fig. 2.31 Construction of Idealized Hydrodynamic Parameter Response Spectrum for a Displacement 
Parameter 
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Fig. 2.32 Construction of Idealized Hydrodynami c Parameter Response 
Spectrum for a Velocity Parameter 
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Fig. 2.33 Elasto-Plastic Resistance-Deformation Relationship 
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Fig. 2.34 Inelastic (Elasto-Plastic) Response Spectra, El Centro Earthquake, 

6 s = - 02 * Wg n,ax /K = °- 00 
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Fig. 2.35 Inelastic (Elasto-Plastic) Response Spectra, El Centro Earthquake, 

8 = .02, K n x /K = 0.006 
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Fig. 2.36 Inelastic (Elasto-Plastic) Response Spectra, El Centro Earthquake, 
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Fig. 2.40 Inelastic (Elasto-Plastic) Response Spectra, El Centro Earthquake, 6 = .02,y = 4.0 
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Fig. 3.1 Comparison of Linear Drag Solution (Linearized by Ground Velocity) and 
Nonlinear Drag Solution. System 1L. El Centro Earthquake, 1940, 

N-S Component. 
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Comparison of Linear Drag Solution (Linearized by Ground Velocity) and Non- 
linear Drag Solution. System 1M. El Centro Earthquake, 1940, N-S Component. 
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Fig. 3.8 Comparison of Linear Drag Solution (Linearized by Ground Velocity) and Non- 
linear Drag Solution. System 4H. El Centro Earthquake, 1940, N-S Component. 
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Fig. 3.9 Comparison of Linear Drag Solution (Linearized by Ground Velocity) and Non- 
linear Drag Solution. System 5H. El Centro Earthquake, 1940, N-S Component. 



206 



♦ 

• 4 



• 4 






• ♦ 
• 4 





LO 


o 


o 


• 


• 


o 


o 


II 


II 


CO 


CO 


CQ 


c n 


• 


4 



>> 

o 

c 

<u 

3 

cr <u 

CL) CD 
S- c 
Ll. fO 
QC 
-£Z 
CD 



CD 

I 

LO 

LO 

I 

I 

CO 

CO 

I 

C\J 

C\J 

I 



>> 

u 

c 

<1) 

3 

cr aj 
0 > CD 

c. c 

Ll. 03 

-o 



CO 

CO 

03 



CD 

LO 

LO 

I 

I 

CO 

CO 

I 

CNJ 

OJ 

I 



co 

CO 

03 



♦ • 



-♦* 



>> 

u 

c 

<D 

3 

cr cu 
<D CD 
S- C 
Ll. 03 
QC 
3S 

O 



CD 

I 

LO 

LO 

I 

^r 

^r 

i 

co 

CO 

I 

C\J 

C\J 

I 



CO 

co 

03 



O 

C\J 



o 

o 



o 

CD 



O 

CO 



(6euo ueauLLuoM) uoliuoislq 6inads # xew 
(6ea'a ueauny uoliuoislq 6uuds ’xew 



Fig. 3.10 Comparison of Linear Drag Solution (Linearized by Ground Velocity) and Non- 
linear Drag Solution. System 6M. El Centro Earthquake, 1940, N-S Component. 
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Fig. 3.11 Comparison of Linear Drag Solution (Linearized by Ground Velocity) and Non- 
linear Drag Solution. System 6H. El Centro Earthquake, 1940, N-S Component. 
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Fig. 3.12 Comparison of Linear Drag Solution (Linearized by Ground Velocity) and Non- 
linear Drag Solution. System 7M. El Centro Earthquake, 1940, N-S Component. 
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Fig. 3.13 Comparison of Linear Drag Solution (Linearized by Ground Velocity) and Non- 
linear Drag Solution. System 7H. El Centro Earthquake, 1940, N-S Component. 
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Fig. 4.1 Response of Towers A, B, C, and D to El Centro 
Earthquake, 1940, N-S Component. 3 = 0.01 
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Fig. 4.1 (Continued) 
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Fig. 4.2 Idealized Elastic Deformation Response Spectra for the First Five 
Modes of Tower C. e = .01 in Each Mode. 
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Fig. 4.3 



Response of Towers A, B, C, and D 
Using Idealized Spectra 
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Fig. 4.3 (Continued) 
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